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Paper 1193 


Journal of the 
PIPELINE DIVISION 


Proceedings of the American Society of Civil Engineers 


DESIGN AND STRENGTH OF 
WELDED PIPE LINE BRANCH CONNECTIONS2 


E. C. Rodabaugh! and H. H. George2 
(Proc. Paper 1193) 


INTRODUCTION 


When an opening is cut into a pipe the uniform distribution of stress under 
pressure is interrupted; the resulting stresses around the opening are signifi- 
cantly higher than those in unperforated pipe. Accordingly, when this opening 
is made for a branch connection, reinforcement is required to obtain strength 
equal to that of the unperforated pipe and to provide for external loadings 
which may be imposed through the branch pipe as shown in Figure 1. 

As yet there is no theoretical analysis of the strength of pipe line branch 
connections, the problem being extremely complex due not only to the shape 
of the structure but also due to the presence of four independent variables— 
the diameters and thicknesses of the run and branch pipes. Consideration of 
laterals adds another variable, the lateral angle. Accordingly, the design of 
branch connections is based entirely on empirical rules derived from experi- 
ence and tests. The empirical rules generally used for the reinforcement of 
fabricated branch connections in pipe lines are given in the American Stand- 
ard Code for Pressure Piping, ASA B31.1. Similar rules for reinforcement 
of branch openings are also given in the ASME Boiler and Pressure Vessel 
Code. 

Brief discussions of loadings, Code requirements and field failures of 
branch connections are presented. A rather extensive amount of test data has 
been accumulated and is reviewed and correlated with Code requirements and 
field experience in the following order. 


1. Pressure loading 
a) Static pressure 
b) Cyclic pressure 


Note: Discussion open until August 1, 1957. Paper 1193 is part of the copyrighted 
Journal of the Pipelines Division of the American Society of Civil Engineers, Vol. 


83, No. PL 1, March, 1957. 
a. Presented at a meeting of the Committee on Pipelines, Construction Divi- 


sion, ASCE, October 16, 1956, Pittsburgh, Pa. 

1. Senior Research Engr., Tube Turns, Div. of National Cylinder Gas Co., 
Louisville, Ky. 

2. Manager of Research Section, Tube Turns, Div. of National Cylinder Gas 
Co., Louisville, Ky. 
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2. External loadings 
a) Static external loads 
b) Cyclic external loads 


Finally, the results of the test data and field experience are summarized 
in a general way under “Recommended Design Practice.” 


Loadings on Branch Connections 


The types of loadings that may be applied to a branch connection are illus- 
trated in Figure 1. Generally in field service a combination of these loads 


will be applied. 
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Fig. 1. Loadings Applicable to a Pipe Line Branch Connection. 


The internal load, i.e. the pressure, exerts a uniform load normal to the 
pipe walls; this load may be negative in the case of vacuum service. The end 
thrusts due to the pressure may or may not be present, depending upon how 
the pipes are anchored adjacent to the intersection. The internal pressure 
may vary due to such factors as: 


a) change in flow rate; 
b) intermittent service; 
c) pulsations set up by the pump and/or line resonance. 


The external loads have been classified as bending, twisting and axial; 
bending being further subdivided into “in-plane” and “out-of-plane.” In ser- 
vice, the bending load may be between these planes; however, for analytical 
purposes, the actual bending load may be separated into its in-plane and out- 
of-plane components. Static external loads are imposed by unsupported 
weight of the piping, installation misalignment, foundation or support settle- 
ment and similar factors. Cyclic external loads may be present due to 
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conditions such as: 


a) cyclic temperature changes producing variations in the pipe line length 
and corresponding variations in loads on the branch connection; 


b) mechanical vibrations of the branch line. 


Code Requirements 


Rules for reinforcement of fabricated branch connections, primarily from 
the standpoint of providing adequate internal pressure strength, are given by 
four sources: 


1. ASA Code for Pressure Piping, Section 6, “Fabrication Details” 


. ASA Code for Pressure Piping, Section 8, Paragraph 831.4, “Reinforce- 
ment of Welded Branch Connections” 


. ASME Unfired Pressure Vessel Code, Paragraphs UG-36 to UG-46, 
“Openings and Reinforcements” 


. ASME Power Boiler Code, Paragraph P-268, “Openings and Reinforce- 
ments” 


The basic concept of reinforcement of branch connections is illustrated by 
Figure 2. Essentially, these rules require that the metal in the run-pipe cut 
out by the opening must be fully replaced by extra metal placed within a zone 
close to the opening. Details of reinforcement requirements and restrictions 
on the opening size, ratio of opening to header size, angle between header and 
branch, etc. vary with the specific Code. All of the Codes point out that the 
rules provide reinforcement for internal pressure; external loadings are not 
evaluated or are assumed to be those normally applied to such connections. 
Cautionary notes are given requiring special consideration where the external 
loads are unusually high or may be cyclic. Special care is also required 
where the pressure may be cyclic. 


“Area of Reinforcement’ Enclosed by —— — — —— Lines. 


Reinforcement Area Required Ag =(4)(t) 


Area Available as Reinforcement = Ay + Ag + A; 


Fig. 2. Typical Code Reinforcing Requirements. 


Branch connections which are sold as “Steel Butt-Welding Fittings” are 
covered by American Standard ASA B16.9. In this standard, adequate strength 
of the fittings is insured by means of proof tests; the standard requires that 
actual bursting pressure of the fittings must be at least equal to the computed 
bursting pressure of the seamless pipe of the size, schedule number (wall 
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thickness) and material designated by the marking on the fitting. Reputable 
manufacturers of such steel butt-welding fittings run numerous hydrostatic 
tests on their pipe fittings to insure that their design is such as to meet this 


requirement. 


Field Failures and Their Causes 


Several failures of saddle-reinforced branch connections on gas and oil 
transmission lines have been reported during the past few years; a failure of 
this type is shown in Figure 3a. 

Another type of commonly encountered field failure involves a small 
branch connection to a relatively large run pipe, reinforced by pads or sad- 
dles or, in some cases, unreinforced; mechanical vibration or movement of 
the branch pipe with respect to the run has caused failures of such connec- 
tions; a failure of this type is shown in Figure 3b. This type of branch con- 
nection formerly gave considerable trouble in steam power plants, however, 
experienced pipe designers in that field now recognize the inherent high 
stress conditions in such connections and make the small branch pipes of 
XXS or Sch 160 wall thickness and, in some instances, maintain a certain 
minimum diameter for small branch connections. 


CRACK AT 
TOE OF 


Fig. 3a Fig. 3b 
Fig. 3. Field Failures. 


The field failures generally involved connections which were adequately 
reinforced according to Code rules. Further, hydrostatic tests on similar 
branch connections would indicate that the connections were adequately de- 
signed for the operating pressure. Since field failures have occurred at 
normal operating pressures, it appears to the authors that such field failures 
of branch connections are not usually due to static pressure load but are 
probably due to such loads as: 


a) cyclic internal pressure 
b) excessive static external loads 
c) cyclic external loads 
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In considering the possible causes of field failures, other tactors besides 
the loads must be considered. A fabricated branch connection involves a con- 
siderable amount of manual welding which introduces the welding quality as a 
factor. The weld around the intersection between the run and branch pipes is 
rather difficult to make with consistently good quality, even for an experi- 
enced pipe welder. The ends of the branch pipe and/or the edge of the hole 
in the run pipe must be carefully beveled if a full penetration weld is to be 
obtained. In welding a pad or saddle type reinforcement, the fillet weld be- 
tween the edge of the pad or saddle and the run pipe can result in an undercut 
at the toe of the fillet weld in the run pipe and, as will be shown later, this is 
a favored location for failure in such reinforced connections. Another factor, 
in connection with welding, is the possibility of starting cracks by welding. 

It is known that, in the case of steels rather high in carbon and manganese 
welded at low ambient temperatures in a constrained position, welding may 
start cracks either in the weld metal itself or in the adjacent base material. 
While the authors are not aware of any field failures directly attributable to 
weld associated cracking, it cannot be entirely ruled out as a possible con- 
tributor to field failures. 

Other factors which may influence or cause field failures are corrosion 
and, at high temperatures, creep or graphitization. 


Static Pressure Tests 


Static pressure tests may be briefly described as tests wherein the pres- 
sure is gradually increased until the branch connection fails by yielding and 
rupture. In conjunction with these static tests, additional information is 
gained by measuring the strains at critical locations as functions of the pres- 
sure. These strains, which may be measured by electrical resistance or 
mechanical strain gages or by brittle coating, may then be converted to 
stresses and compared with allowable stresses for the material under static 
or cyclic loading conditions. Either strains or overall dimensional changes 
may be plotted against the applied pressure to determine the pressure cor- 
responding to yielding; however, in the case of a complex shape such as a 
branch connection, the practical significance of yielding has not been 
evaluated. 

It is impossible here to describe and give the results of all of the static 
pressure tests that have been run on branch connections. A list of references 
to published test data is given in the Appendix, along with a brief abstract of 
their contents. In the following we will cover typical tests on the more widely 
used types of branch connections and indicate in general terms the signifi- 
cance of such tests. 


a) Unreinforced Branch Connections 


For the present purposes, an unreinforced branch connection will be de- 
fined as a branch connection made by welding a branch pipe to the opening in 
the run pipe without the provision of additional reinforcement. It is possible 
that such a connection may be fully reinforced according to Code rules. This 
occurs where the branch pipe is very thick in relation to its diameter; also, 
if the run pipe is twice as thick as is required to hold the design pressure, 
the branch connection is generally fully reinforced for the design pressure 
according to Code rules. 
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Figures 4a and 4b show failures produced by hydrostatic tests to rupture 
i.e. burst tests. The left hand photographs show failures in the “crotch” of 
the intersection, which is the more normal location of the rupture; however, 
rupture often occurs at the sides as shown by the right hand photographs; the 
intersection weld being a weak point. Table 1 describes and gives the results 
of a number of typical burst tests of unreinforced branch connections. 


Table 1: Typical Test Results, Burst Tests of Unreinforced Branch Connections 


+50 = 
610 


2375 


* P is tne actual ourst pressure of an equivalent straight 
pipe, when availa.le, otherwise it is the calculated burs’ 
pressure of equivalent pipe usin: the tarlow formla with 
run pipe dimensions and ultimate tensile streagth of run 
pipe material. 

® avera;e of three tests. 


Based on hydrostatic burst tests of unreinforced tees and laterals, an em- 
pirical formula for the strength of such branch connections has been 
developed which is: 


1.5 


(1-0.7E sin ) (1) 


where P, = burst strength of branch connection, psi 
Pc = burst strength of uncut pipe of same O.D. and wall thickness as 
the run of the branch connection, psi 


d = branch inside diameter, in. 

dr = run inside diameter, in. 

E = weld efficiency; E = .8 for field welding; E = 1.0 for shop weld- 
ing where procedures for insuring good quality welds are pro- 
vided. 

@ = angle between run and branch (acute angle for laterals) 


In addition to hydrostatic tests to rupture, a number of tests have been 
made on unreinforced branch connections in which stresses around the branch 
were determined by means of strain gages. As might be expected, stresses in 
an unreinforced branch connection are highest at the junction between the 
branch and run pipe, the stresses rapidly damping out with distance from the 
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intersection. Variation of stresses around the intersection of three unrein- 
forced branch connections are shown in Figure 5, abstracted from Reference 3. 
The ratios of maximum measured stresses to the nominal calculated stresses 
are plotted against the position around the intersection. 

In the absence of a theoretical analysis, it is difficult to draw any general 
deductions as to stresses in unreinforced branch connections. It should be 
noted that stresses determined by strain gage measurements are valid only 
as long as the structure behaves elastically. In particular the stresses so ob- 
tained do not necessarily give any information as to the burst pressure of the 
connection. For example, the strain gage measurements of the 24" x 24" 
fabricated unreinforced branch connection indicated a maximum stress which 
was about five times the stress in the run pipe at locations remote from the 
branch connection; however, this particular branch connection, when sub- 
jected to a burst test, had about 75% of the strength of the run pipe or a 
“stress ratio,” at the burst pressure, ot only 1.33 instead of 5. This is due 
to the fact that when a branch connection is fabricated from a ductile material, 
the material can yield and readjust to a more favorable shape before breaking. 
On the other hand, when loads are applied cyclically the stresses indicated by 
strain gages are significant. As will be discussed later herein, a rough cor- 
relation exists between static strain gage measurements with the number of 
cycles to failure and failure location in cyclic pressure tests. 


b) Fabricated Reinforced Branch Connections 


A fabricated reinforced branch connection is defined herein as a branch 
connection which has additional reinforcement components placed around the 
branch outlet; these reinforcement components usually consist of saddles, 


pads, encirclement sleeves or combinations thereof. Table 2 describes and 
gives the results of a number of typical burst tests of reinforced branch con- 
nections. Typical failures are shown in Figures 4c, 4d and 4e. 


Table 2: Typical Test Results, Burst Tests of Reinforced Branch Connections 


% Reinforcing Burst Failure Tests 


+«éj4 

Branch per ASA Pressure, psi 
rs 

7150 
3750 
6350 
2350 
2100 


4.600 


8750 


1865 


2490 
1440 
2060 


* 411 tests were on 90° branch connections (tees). 
® See figures 7 and § for sketches of these types of reinforcements. 


> Percent of reinforcing required to make intersection equal in strength to the run pipe in accordance with 


the American Standard Code for Pressure Piping, ASA B31.1-1955, Par. 639. 


* Failures of pad or saddle reinforced branch connections usually occur on the side at the toe of the fillet 


weld between reinforcement and run pipe as shown in Fig. kc. 
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Adequately reinforced branch connections, such as those tested, fail in the 
pipe at about the pressure at which unperforated pipe would be expected to 
fail. In other words, tests available to date indicate that a 90° branch con- 
nection adequately reinforced in accordance with ASA Code for Pressure Pip- 
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Fig. 5. Experimentally Determined Stresses in an Unreinforced Branch 
Connection (Abraham and McClure). 
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ing rules is equivalent to the run pipe in burst strength. It will be shown 
later; however, that these types of reinforced branch connections can be much 
weaker than straight pipe if the pressure is cyclically applied. 


c) Manifold Welding Fittings 


Manifold welding fittings are used where several branch outlets in a line 
are needed in close proximity to each other. A photograph of a typical mani- 
fold welding fitting is shown in Figure 6. This type of fitting is made by draw- 
ing one or more branch outlets in a section of run or header pipe, the header 
pipe being considerably thicker than required for unperforated pipe; in this 
design the reinforcement is integral with the fitting. Welding of the branch to 
the run is accomplished by means of a relatively simple circumferential butt 
weld thus avoiding the welding difficulties inherent in the complex branch-to- 
run pipe weld encountered in fabricated branch connections. 

The failure resulting from a burst test of a section of a typical manifold 
fitting is shown in Figure 4f; in this particular case the branch pipe was the 
weakest member. This manifold fitting, a 10" x 6" size, was designed to be 
equal to burst strength to 10" x .500"" wall, Al06 Gr B run pipe. The burst 
pressure was 5810 psi which is 1.41 times the calculated burst pressure of 
the equivalent pipe, thereby proving adequacy of the reinforcement in ac- 
cordance with ASA B16.9 requirements. 


d) Butt-Welding Forged Tees 


Butt-welding tees manufactured by the authors’ company are forged in a 
manner causing a preferential increase in wall thickness at the locations of 
highest stress. In addition, the “barrel” shape of the run and the large crotch 
radii combine to minimize stress concentrations from either pressure or ex- 
ternal bending loads. The failure resulting from a burst test of a typical 
forged tee is shown in Figure 4g. Since the tee is designed to be stronger 
than the matching pipe, failure occurred in the run pipe in this particular 
case at a pressure of 2600 psi which is 1.19 times the calculated burst pres- 
sure of the 18" O.D. x .375" wall, A106 Grade B matching pipe. In order to 
obtain information on the burst strength of the forged tee itself, it is usually 
necessary to either use heavier pipe legs or to reinforce the pipe legs with 
steel bands as shown in Figure 4h. In this case failure occurred at 5950 psi 
which is 1.34 times the calculated burst pressure of the 6.625'" O.D. x .280" 
wall, A106 Grade B matching pipe. 


Cyclic Pressure Tests! 


As the result of several failures of saddle reinforced branch connections 
reported in gas and oil transmission lines, serious thought was directed to- 
wards improving the design of the reinforcement of such openings. Since 
failures usually occurred at the toe of the weld between the saddle and run 
pipe (see Figure 3a), attention was directed towards a type of reinforcement 
which would completely encircle the run pipe. Several varieties of these en- 
circlement type reinforcements were proposed; these are shown in Section 8 
of the ASA Code for Pressure Piping, ASA B31.1, 1955. Some of these rein- 
forcements were designed for application to “hot tap” connections, wherein a 


1. Part of the test work described in this section has been previously pub- 
lished in Reference (7). 
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branch connection is installed without stopping flow through the main pipe 
line; they may also be used for normally constructed branch connections. To 
assist pipeline operators in evaluating the respective merits of these new en- 
circlement type reinforcements, and to generally evaluate the strength of 
branch connections, the authors’ company initiated an extensive test program 
in which the relative strength of the various branch connections was deter- 
mined by application of cyclic pressure. 

General views and enlarged sections of significant details of the test 
specimens are given in Figures 7, 8 and 9. Figure 7 shows the first series 
of test specimens which all involve 22" O.D. x 5/16" wall run pipe and 
10.750"' O.D. x .500"' wall branch pipes. Figure 8 shows the second series 
of test specimens which involve 24" run pipe with 4", 8"' or 12" branches 
with various reinforcements.2 Figure 9 shows the third series of test speci- 
mens which are sections of manifold welding fittings. In the first series ap- 
proximately half-scale models of the variants shown were constructed and 
tested in addition to tests on full scale specimens. The material used in the 
models was different from the full scale specimens only in that the run pipe 
was A106 Grade B pipe instead of API-5LX-52 pipe. 

The test apparatus developed for the cyclic pressure tests is shown in 
Figure 10. Briefly, the operation was as follows. The specimen was water 
filled and air-vented and then connected through a header to the discharge 
end of a reciprocating pump. The pump was started; as soon as the pressure 
reached a pre-selected maximum, a pressure control switch opened a by- 
pass valve allowing discharge into a storage tank connected to the pump suc- 
tion; as soon as the pressure dropped to a pre-selected minimum, another 
pressure control switch closed the by-pass valve allowing the pressure to 
rise again and the cycle was repeated. Each cycle operated a counter. 
Failure of a specimen caused the pressure to drop below the pre-selected 
minimum, at which time a low pressure switch automatically shut down the 
entire apparatus, the counter reading indicating the cycles to failure. 

In general, the upper pressure limit was set around 90% of the calculated 
yield pressure of the run pipe; lines are often tested to this pressure and it 
might also be attained as a result of pressure pulsations. The lower pres- 
sure used was of the order of one-half the yield pressure. While pressure 
variations of this order of magnitude are probably not frequent in actual ser- 
vice, they were considered necessary to produce failure in a reasonable 
length of time. This experimental technique of applying intensified loading 
is, of course, widely used and accepted as a basis of obtaining data on com- 
parative performance. 

The tests were continued until failure of the test specimen occurred or 
until a multiple of the cycles causing failure of a corresponding saddle rein- 
forced connection had been applied. 

A first quantitative yardstick of the relative performance of the different 
variants tested is given by cycles each sustained before failure; these are 
shown on Figure 12. Since several sizes were involved, it is perhaps more 
informative to relate the cyclic life of the individual variants to that of a sad- 


2. Series Il specimens were used at Battelle Memorial Institute for compre- 
hensive strain gage tests prior to being subjected to cyclic pressure tests 
by the authors’ Company. This was part of an extensive investigation of 
stresses in pipeline branch connections sponsored by the American Gas 
Association. (8) 
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Fig. 7. Specimens for Series I Cyclic Pressure Tests. 
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Fig. 9. Specimens for Series III Cyclic Pressure Tesis. 


Fig. 10. Cyclic Pressure Test Apparatus. 
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dle reinforced branch connection. The saddle reinforced connection was 
chosen as a basis for evaluation since it was the type which had occasionally 
failed hence this type reinforcement may be considered as a border-line 
design. Branch connections which are significantly superior to the saddle 
reinforced design would presumably be adequate under field conditions; those 
connections which are significantly inferior to the saddle reinforced design 
might be expected to fail under severe field service conditions. 

In order to generalize the available test data on saddle reinforced branch 
connections, a plot of equivalent stress versus cycles to failure was con- 
structed as shown in Figure 11. The cyclic pressure tests may be con- 
sidered as applying a mean stress upon which a cyclic stress is superim- 
posed. The relation between completely reversed cycles of stress and a 
combination of mean stress plus cyclic stress was adopted from the relation 
suggested by Moore and Kommers.( 


S,=S_+1/38, (2) 


where S, = equivalent stress for completely reversed cycles 
Sy = variable stress component 
Sc = constant stress component 


Following the type of analysis used for cyclic bending tests,! it was found 

that the relation between stress and cycles-to-failure could be represented 
by an equation of the type Se = CN--2; for the particular proportions of sad- 
dle reinforced branch connections tested, a value of C = 216,000 fits the test 


data fairly well. 
SIZE 
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Fig. 11. Saddle Reinforced Branch Connections, Stress vs Cycles to Failure. 
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Of the seven saddle reinforced connections tested; five are roughly geo- 
metrically similar. The other two (24" x 4" and 24" x 8", specimens C5 and 
Cg), since they involve a smaller ratio of branch to run size, would not, a 
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Fig. 12. Comparison of Cyclic Lives of Variants Tested. 
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Fig. 13. Failure of Rolled-and-Welded Pipe. 
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priori, be expected to fall on the same S-N curve as the other five speci- 
mens. It was found, however, that the one line on Figure 11 fairly well rep- 
resents the data for all seven specimens. Apparently the controlling factor 
is the discontinuity at the edge of the saddle-to-run pipe weld, since failure 
consistently occurred there, and possibly the geometry of the rest of the con- 
struction, at least over a limited range, has little influence over the control- 
ling stress. This hypothesis is somewhat supported by the pad reinforced 
connection tests, which failed at the same location at about the same number 
of cycles. However, the generalizations are based on rather limited data and 
extrapolations to shapes or proportions radically different than those tested 
may not be justified. 

Having established a reference line for saddles, the relative life of various 
types of branch connections are then as shown in Figure 12, wherein the dif- 
ferent types of branch connections are arranged in order of decreasing life. 

As anticipated, the straight pipe (A) lends itself well as an upper aim point 
and the unreinforced intersection (B) sets a lower performance limit. The 
ordinary saddle reinforced connection (C), our yardstick, because it is the 
construction all other designs are intended to improve upon, lies in between. 
The forged tee (K) did not fail in a single instance and thus justifies its 
recognized position as the full equivalent of the pipe it is intended to match. 
The manifold welding fitting (L), encirclement saddle welded along the sides 
(G) and the encirclement sleeve-saddle (E) all had cyclic lives essentially 
equal to that of rolled-and welded straight pipe. The bolted encirclement 
saddle (H) is an improvement upon the ordinary saddle reinforcement but not 
the equivalent of straight pipe. The encirclement saddle welded at the crotch 
and back (F) was somewhat erratic and, on the average, behaves only slightly 
better than the ordinary saddle reinforcement. The pad reinforced connec- 
tions (M) and the encirclement sleeve (D) gave about the same performance 
as the ordinary saddle reinforced connections. Finally, the encirclement tee 
(J) has been found erratic in performance and generally inferior to the 
ordinary saddle, contrary to expectations. 

Failures of the rolled and welded steel pipe consisted of cracks in the 
longitudinal direction in the vicinity of the longitudinal seam weld in the pipe. 
Generally the split paralleled the seam weld, following one of several shallow 
grooves observable in the outer surfaces of the pipe. In one failure where the 
split was short, it was noted that both the length and width of the crack was 
greater on the outside than on the inside clearly indicating that the failure 
started on the outside. Figure 13 shows an overall view and close-up of such 
a typical failure in resistance weld pipe. The large spread in cyclic life be- 
tween straight pipe specimens, one failing at 144,000 cycles and another run- 
ning for 900,000 cycles without failure, probably reflects the variation in the 
severity of the mechanical or metallurgical notches. 

The unreinforced connections failed either by cracking across the weld in 
the crotch, or along it at the side, both types of failure being experienced 
Simultaneously in one assembly; see Figure 14. This is in accord with strain 
gage test results obtained by other investigators; for connections where the 
branch size is one-half that of the run or less, the maximum stresses are 
nearly uniform around the perimeter of the intersection, only the direction 
changes. 

In case of the connection reinforced with an ordinary saddle, the run pipe 
invariably split along the side at the edge of the saddle as shown in Figure 15. 
A cross section showing the start of a crack on the side opposite to that where 
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14a. Arrows Show 14b. Detail of Crack 14c. Detail of 
Location of Cracks in Crotch Crack in Side 


Fig. 14. Failure of Unreinforced Connection. 


the rupture occurred is illustrated by Figure 16; this section reveals how a 
fatigue crack starts at the location of marked chance in section at the toe of 
the attachment weld, which discontinuity in some instances may be accentuat- 
ed by under-cutting during welding. Noticeable change in texture of the 
cracked face evident in the view in Figure 15 indicates how far the fatigue 
crack progressed inwardly through the wall before this was weakened so that 
it could no longer sustain the internal pressure and burst open. 

The pad reinforced connections failed in the same manner as the saddle 
reinforced connections as shown by the typical failure in 17. 

The encirclement sleeve reinforced connection failed by splitting through 
the crotch as shown in Figure 18. Both the test variant itself and the failure 
location bear close resemblance to the unreinforced intersection. By the ap- 
proximate correlation between stress and cycles to failure suggested by 
equation (4), a comparison of cycles to failure of the sleeve-reinforced con- 
nections with the unreinforced connections implies that the 1/2"' thick sleeve 
in the 22" x 10"' assembly served to reduce the controlling stress to 55% of 
the value without this reinforcement; the 3/8" thick sleeve on the 24" x 12" 
assembly similarly reduced the stress to 75% of the value without reinforce- 
ment. These would appear to be reasonably estimates of the reinforcing ef- 
fect of the sleeves. 

The sleeve-and-saddle reinforcement gave excellent results; two speci- 
mens failed in the pipe run outside of the reinforcement and the third in the 
pipe run under the sleeve, the fatigue crack following the edge of the weld line 
in the run pipe. 

The three variants of the encirclement saddle are very similar to each 
other in intersection detail. Of these three, only the design wherein the en- 
circlement halves were welded along the sides (G) consistently gave results 
in line with those obtained for straight pipe. The performance of Variant F, 
which is distinguished from Variant G by having the encirclement halves 
joined together by lap weld at the back and a butt weld through the crotches, 
was somewhat erratic. In the first model test a leak started under the rein- 
forcement at 47,200 cycles; in order to be able to continue the test, both ends 
of the reinforcement were fillet welded to the run pipe but the assembly only 
sustained an additional 5,500 cycles before the crotch weld joining the two 
halves cracked. This somewhat difficult weld had been made in the normal 
way, without backing. To insure a sound weld, a thin curved backing strip 
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Fig. 15. Failure of Fig. 16. Saddle Reinforced 
Saddle Reinforced Connection, Side Opposite Failure 
Connection. Showing Start of Crack. 


Fig. 17. Failure of Fig. 18. Failure of Sleeve 
Pad Reinforced Connection. Reinforced Connection. 


was used at this location in constructing the remaining three specimens 
of this variant. Probably as a result of this provision, the second model 
4 ran 264,000 cycles without failure; one of the full scale specimens also 

; showed a fair life, while the other performed only moderately better than 
a the ordinary saddle construction, both failing in the intersection weld 
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between the run and branch pipe in the 
crotch area under the reinforcement. 
A typical failure is shown in Figure 19. 
Possibly the reason for the difference 
observed in the performance of vari- 
ants G and F might be explained by the 
hypothesis that shrinkage of the two 
side welds in Variant G in some way 
produces a better fit-up in the critical 
zone than is the case where the welds 
are made at the back and crotch as in 
Variant F, however there is no direct 
evidence that this is true. 

The bolted variant (H) provided 
moderate but consistent improvement 
over the ordinary saddle construction; 
as in the other types of encirclement 
reinforcements failure occurred in the 
intersection weld in the crotch area, 
see Figure 19. With the thought that 
close contact of the saddle with the 
branch-to-run intersection weld might 
have a beneficial effect, intentional variables in the detailed fit-up in this 
location were introduced but no marked effect on relative life of the different 
specimens was produced. This design evolved from the thought that bolts 
should provide an excellent means of bringing the reinforcement into force- 
able contact close to the critical intersection weld; in addition, there is no 
welding against the run pipe such as that which occurs where the sleeve 
halves are joined by a butt weld. In the design actually tested however, the 
bolts were placed too far from the shell and the bolting flanges were some- 
what too light. Improvements in detail of this design might bring its perform- 
ance up so as to be comparable with that of straight pipe. 

The encirclement tee is basically different than the other reinforcements 
previously discussed in that the reinforcing member directly carries the in- 
ternal pressure; the branch pipe is butt welded to the reinforcement and is not 
welded to the run pipe. This type of branch connection gave unexpectedly poor 
results. With the exception of the initial failure in the first model, where the 
butt weld joining the two halves cracked in only 6500 cycles, all failures 
occurred in the fillet welds between the tee and run pipe, the location in all 
except one instance being directly in line with the crotch. As can be seen in 
Figure 20d, the cracks ran in a circumferential direction and apparently 
started at the root of the fillet weld. 

The encirclement tee assemblies, including the fillet welds in all but the 
last specimen, had been proportioned to parallel the assembly used by Del 
Buono, Vissat and Williams(9) in their exhaustive and evidently very carefully 
conducted strain gage tests. It is noteworthy that failure in the cyclic pres- 
sure tests occurred in the precise location where these investigators had ob- 
tained the highest strains, and that the crack closely resembled that found by 
them upon sectioning through the crotch plane after bursting their test 
assembly. 

Since the fillet weld quite obviously constituted a point of weakness, the 
thought naturally presented itself that the performance of this design could be 


Fig. 19. Failure of Encirclement 
Saddle Reinforced Connection, 
Quarter Sectioned to Show Failure 
on Inside of Crotch. 
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20a. General View 20b. Close-up 20c. Same as 20b After 20d. Cross-Section 
Before Grinding Grinding Weld, Crack Through Failure 
Off Weld . Brought Out by Area 
Penetrant Dye 


Fig. 20. Failure of Encirclement Tee Connection. 


Fig. 21. Failure in Manifold Welding Fitting. 


improved by increasing the weld size. Accordingly a special test assembly 
was made by using what was considered a maximum practical weld size. The 
weld leg dimension was increased by more than 50% and the tee end was 
machined with a reverse bevel of 600 and a 3/16" root radius, all of which 
resulted in over 150% increase in the amount of weld metal deposited. This 
design with a heavier weld (Specimen J5) lasted 3 1/4 times as long as one 
with the original weld proportions tested under otherwise identical conditions, 
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Fig. 23. View of End Bank of Fatigue Test Machine Showing a Fabricated 
Tee Mounted for In-Plane Bending Test. 
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however this specimen still had a cyclic life only slightly better than a cor- 
responding saddle reinforced connection. 

Finally we come to the forged tee and manifold welding fittings, these be- 
ing basically different from the preceding variants in that they cannot be used 
for hot tap connection. In none of the three specimens of the forged tee did 
failure occur in the tee itself. The two model tests were run to a large mul- 
tiple of the cycles causing failure in a corresponding ordinary saddle rein- 
forced connection; in the full scale specimen the test was terminated by 
failure of the run pipe with the tee left unharmed. It is apparent that the in- 
tegral reinforcement, large crotch radii, preferential thickening in highly 
stressed areas and gradual change in wall thickness make the forged tee the 
full equivalent of the pipe it is designed to match. 

The manifold welding fittings also gave a good account of themselves but 
they were not equal to the forged tees. All three manifold welding fittings 
had a cyclic life of the same order as that obtained for straight pipe. The 
failures consisted of cracks in the crotch area; a typical failure is shown in 
Figure 21. It was not expected that the manifold welding fittings would be as 
good as the forged tees since they do not have the generous crotch radii and 
preferential thickening in highly stress areas which are obtained with forged 
tees. 

Results of the cyclic pressure tests are correlatable with strain gage 
tests. In the Series II tests, wherein the specimens had been used for strain 
gage tests. In the Series I tests, wherein the specimens had been used for 
strain gage tests prior to the cyclic pressure tests, all seven specimens 
failed in locations at or near the points of maximum measured stress. As 
was pointed out earlier, failures of the encirclement tees (J) occurred at the 
location of maximum measured stress. A rough correlation between stresses 
and cycles to failure was also observed; since strain gages do not measure 
the effect of local notches, such as weld defects or undercutting, it would not 
be expected that good correlation between stress measured by strain gages 
and cycles to failure would necessarily exist. 

Although all of the reinforced fabricated branch connections used in the 
cyclic pressure tests were more than 100% reinforced by the ASA Code for 
Pressure Piping Rules, these connections were not all equal to the longitudi- 
nally welded straight pipe in cyclic life. Accordingly, while it was concluded 
that 90° fabricated branch connections reinforced per the ASA rules were 
equal to straight pipe in burst strength, the same conclusion cannot be drawn 
for cyclically applied pressure. 


Static External Load 


Figure 22 shows a test in which a combination of static internal pressure 
and static out-of-plane bending load was applied to produce rupture. In this 
test, a 22" x 10" saddle reinforced branch connection, constructed as shown 
in Figure 7, was pressurized with water to 1200 psi; then an out-of-plane 
bending force was applied and gradually increased until failure occurred at 
a bending moment of 300,000 ft-lb. The pressure corresponded to a circum- 
ferential stress in the run pipe of approximately 72% of its yield strength; 
the bending moment at failure gives a calculated maximum stress in the 
branch pipe of 73,000 psi. This test indicates that very large static bending 
loads can be carried by a reinforced branch connection. 

Stresses resulting from static bending moments are shown in Figure 5. 
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Cyclic External Loads 


Cyclic bending tests, in which about one hundred different branch connec- 
tions were subjected to cyclic bending loads to failure, were run using the 
test machine shown in Figure 23. In these tests, one end of the assembly was 
anchored rigidly and the other, a hinged end, was forced to deflect cyclically 
through a given amplitude maintained by the machine eccentric setting. The 
specimens were subjected to individual load-deflection calibrations, filled with 
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Fig. 24. Reinforced Fabricated Branch Connections Used in Cyclic Bend- 
ing Tests. 


water to provide a ready means of detecting failure, and finally flexed cycli- 
cally through a predetermined amplitude until a leak signalled a crack through 
the entire wall. This was done at several amplitudes for each type of speci- 
men so that upon completion of tests a stress vs cycles-to-failure plot (S-N 
curve) could be constructed. The stress is the nominal stress in correspond- 
ing size and wall thickness straight pipe calculated by the ordinary beam 
formula: 


(3) 


= nominal bending stress, psi 
applied load (lb), obtained from load-deflection calibration 
L = lever arm from point of load application to point of failure, in. 
Z = section modulus of corresponding straight pipe, cu-in. 


Details of reinforced test specimens are shown in Figure 24; in addition, 
various designs of forged butt-welding tees and unreinforced straight fabri- 
cated branch connections with dimensions as indicated on Fig. 25 were tested. 
Results of these tests as S-N plots are shown in Figure 25, along with results 
of bending fatigue tests by Blair. (1) Photographs of typical failures are shown 
in Figure 26. 

Tests conducted by Blair were run with static internal pressure applied to 
the connections while applying the cyclic bending loads. Accordingly, in these 


The data given herein on cyclic bending tests is principally abstracted from 
reference (5). For more complete details of these tests and their evaluation, 
this paper should be consulted. 
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tests there was a constant stress due to pressure as well as cyclic stresses 
due to bending; for these specimens an arrow pointing upward is shown, the 
tip of the arrow indicating the equivalent stress for complete reversal in ac- 
cordance with equation (2). 

It is pertinent to note the location of failures in the forged tee, indicating 
the areas of highest stress. It is in just this highly stressed area where the 
intersection weld of a fabricated branch connection is located. Since it is 
practically impossible to make such an intersection weld without some de- 
fects and notches, such fabricated connections, even though heavier than the 
seamless forged tee, are liable to premature fatigue failure due to the stress 
concentrations caused by weld irregularities. Pad reinforced intersections 
also failed along the intersection weld. In the case of the saddle reinforced 
connections, the saddle was able to support the intersection weld so that the 
failure location is transferred to the outer edge of the saddle at the toe of the 
fillet weld between saddle and pipe. In this location, unless care is used, 
fillet welding can cause an undercut into the run or branch pipe, thereby 
creating an additional notch in a highly stressed area. The welding of these 
tests specimens was done with great care under ideal shop conditions; with 
carelessly made welds or with poor welding conditions the fatigue life might 
be very adversely affected. 

Similar bending fatigue tests were run on circumferentially butt-welded 
straight pipe and the results of these approximately 50 tests were also plotted 
on an S-N graph. The relative strengths of the various types of branch con- 
nections were evaluated in relation to the strength of such butt-welded pipe. 
It was found that the test data could be reasonably well represented by the 
formula: 


2 


iS = 245,000N ° 


(4) 


where S = nominal applied bending stress in corresponding straight pipe, 
psi 
N = cycles to failure 
i = stress intensification factor; equal to unity for circumferential- 
ly butt-welded straight pipe 


The values of the stress intensification factors in relation to butt-welded 
pipe are shown in Figure 25. In attempting to evaluate and generalize these 
test results, it was noted that failures in forged tees were analogous in 
location and direction to those occurring in elbows or curved pipe submitted 
to the same type of fatigue test. Stresses in curved pipe can be readily 
evaluated by theoretical methods; accordingly the theory of curved pipe was 
used as a basis for generalizing the experimental fatigue life data on branch 
connections, leading to the formulas shown in the table 3. 

The stress intensification factors based on these tests have been incor- 
porated in the American Standard Code for Pressure Piping, along with fac- 
tors derived from similar tests on other piping components such as elbows, 
mitres and flanged joints. 

Each of the individual cyclic bending tests represents considerable ex- 
penditure of labor, hence the number of specimens and variants were neces- 
sarily limited. Accordingly, the generalizations derived are based on rather 
limited data; only the pressing need for design information induced the 
attempt to translate the imperfectly defined pattern of results into working 
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Table 3: Formulas for Stress Intensification Factors 


DESCRIPTION SKETCH h STRESS INTENSI- 
FICATION FACTOR 


Full-size 44 


Welding Tee 
per ASA B16 9 (1) 


r 
t 


x— 


Full-size pod 
or Saddie Reinforced 
Fobricoted Tee (i) (2) 


Full-size 
Unreinforced 
Fabricoted Tee (i) 


() Also use for reducing-outle tees in absence of directly applicable doto. 
2) When Tis greater thon |.5t,use h= 4.05 +. 


constants. As a corollary, extrapolation of the data given to shapes or pro- 
portions radically departing from those tested is liable to appreciable error. 


Recommended Design Practice 


The test data on branch connections accumulated to date represent an im- 
mense amount of time and effort; nevertheless, it is not possible to quantita- 
tively apply the test data to branch connections in general. The test data does, 
however, indicate some general design practices which should be followed. 

The design should first satisfy the rules of the Codes cited herein. If ser- 
vice conditions are more severe than contemplated in Code requirements, i.e., 
if cyclic pressure, high static bending loads or cyclic bending loads of any 
magnitude may be applied, then consideration should be given to the selection 
of a suitable type of branch connection which will not only meet Code require- 
ments but will be of optimum design for the expected service conditions. 

In an optimum design, one must balance the first cost of the installation 
against the possibilities of failure and the consequence of such a failure. In 
some pipe lines which may be shut down with little effect on associated pro- 
cesses, a failure may mean only the necessity of replacing the failed com- 
ponent. In other pipe lines, however, a forced shut down may be extremely 
costly. Failure of lines carrying lethal or inflammable fluids or gases at 
high pressure may result in great economic loss and hazard to life. 

The test data and field experience all tend to indicate that branch connec- 
tions with integral reinforcing, large intersection radii, gradual changes in 
wall thickness and preferential thickening in highly stressed areas will give 
optimum strength. That this should be true of branch connections is not sur- 
prising; parallel experience with machine design has long ago led to the use 
of generous radii and gradual thickness changes to reduce the possibility of 
fatigue failures. Branch connections with such characteristics in varying de- 
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Fig. 26. Failures in Cyclic Bending Tests. 
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grees are commercially available as forged tee and crosses, manifold weld- 
ing fittings and drawn outlet laterals. In welding such fittings into a pipe line, 
care must be taken to obtain a full-penetration weld; however, since the weld- 
ing is done in a uniform circumferential joint, this is relatively easy to do as 
compared to making a complicated intersection weld between two pipes. 
Furthermore, the welds in such integral fittings are not ordinarily located in 
the area of highest stress. 

While integrally reinforced branch connections usually give optimum 
strength, there are occasions where it is either necessary or desirable to 
use a fabricated branch connection; e.g. where the size is such that forged 
fittings are not commercially available, for “hot-tap” connections or where 
it is believed that the installed cost of a fabricated connection will be less 
than that using a commercial fitting and that the service conditions, both from 
the standpoint of severity of loadings and consequences of failure, are mild. 

Where fabricated branch connections are to be used, the following types 
of construction are recommended: 
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a) For severe service conditions—use a complete encirclement type of rein- 
forcement such as one of the types shown to be adequate by test. 


b) For moderate service conditions—use saddles, pads or other reinforcing 
members designed to meet Code requirements. Sharp angle laterals may 
require more reinforcing than implied by Code requirements. 


c) For mild service conditions—Code requirements may be met by “unrein- 
forced connections,” however extreme caution should be used to assure 
that substantial external loadings will not be applied, even though the pres- 

sure is of a low order. 


In constructing fabricated branch connections, great care must be taken 
with welding in order to assure adequate service life. In making the branch- 
to-run pipe intersection weld, the ends of the branch should be properly 
shaped and beveled and/or the opening in the run pipe shaped and beveled so 
that full penetration of the weld can be obtained. Lack of penetration of this 
weld not only produces a severe mechanical notch in a highly stressed loca- 
tion but also leaves a crevice which is ideal for promoting cell type corrosion, 
should the fluid have any tendency towards this type of corrosion with the 
particular metal used in constructing the pipe line. The fillet weld between a 
reinforcing member, such as a saddle or pad, and the run pipe must also be 
carefully made to avoid undercutting into the run pipe at the toe of the fillet 
weld. 
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APPENDIX 


Bibliography and Abstract of Published Research Work 
Applicable to Pipe Line Branch Connections 


The following is a list of published papers which have come to the authors’ 


attention and which are believed to be particularly applicable to the design of 
pipe line branch connections. 


1 


“Effect of Openings in Pressure Vessels” by J. Hall Taylor, and E. O. 
Waters, presented at the June, 1933 meeting of the ASME. Hydrostatic 
tests on pressure vessels with 42" inside diameter LD. x 1" wall thick- 
ness with ten openings of various sizes and types of reinforcements. 
Stresses around the openings was determined by strain gage readings on 
the outside surface. Results are presented in graphical form. 


“Investigations of Stress Conditions in a Full-Size Welded Branch Connec- 
tion” by F. L. Everett and Arthur McCutchan, Transactions of the ASME, 
July, 1938. Hydrostatic tests of two 8” x 8” pipeline branch connections; 
one assembly was unreinforced; the other was reinforced with a pad. 
Strains were measured on the outside surface in the area of the intersec- 
tion using a Huggenberger mechanical strain gage. 


“Bursting Pressure Tests of Welded Pipe Headers Show Need for Nozzle 
Reinforcements” by Eric R. Seabloom, Valve World, July-August, 1941. , 
Hydrostatic tests of 8" x 8" pad reinforced branch connection, 12" x 12 
pad reinforced branch connection and a 12" x 6" unreinforced connection. 
Deformations were measured with indicating dial gages. Pressure was 
increased until the test assembly fractured. 


4. “Design of Fabricated Plate Steel Tees, Laterals and Wyes of Large 


Diameters for the Pressure Aqueduct of the Boston Metropolitan District 
Water Supply Commission” by Chester J. Ginder and Edwin B. Cobb, 
Boston Society of Civil Engineers, Vol. 28, 1941. Describes a method of 
designing wye branch connections and laterals in the size range of 60 to 
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108". Details of the designs of several typical branch connections are 
given. 


5. “Strengthening of Circular Holes in Plates Under Edge Loads” by Leon 
Beskin, Journal of Applied Mechanics, September, 1944. Theoretical 
stress distributions are determined around strengthened circular holes 
in plates submitted to edge loads at infinity. Various proportions of cir- 
cular strengthenings are considered. Theory is approximately applicable 
to branch connections in which the run pipe is very much larger than the 
branch pipe so that locally, the run pipe can be considered as a flat plate. 


6. “Stresses in a Cylindrical Shell Due to Nozzle or Pipe Connection” by G. 
J. Schoessow and L. F. Kooistra, Journal of Applied Mechanics, June, 
1945. Tests on a 56” diameter shell with 11 3/4 “ ©. D. unreinforced 
connection. The shell was made in two halves, one half being 1.30" thick 
the other half 2.08“ thick. The branch pipes were 7/8 " thick. The shell 
was anchored while bending or axial loads were applied to the branch 
pipes. Stresses due to these loadings were measured by means of electri- 
cal resistance strain gages. 


7. “Reinforcement of Branch Pieces” by J. S. Blair, Engineering (London), 
July 5, Sept. 6, Nov. 29, Dec. 6, 13, 20, 27, 1946. Hydrostatic, static 
bending and cyclic bending tests of a large variety of branch connections. 
The tri-form type of reinforcement is extensively investigated and 
methods for the design of this type of reinforcement are presented. In the 
hydrostatic tests, the yield pressure was determined by means of a brittle 
coating following which the pressure to cause fracture was determined. 


8. “Experimental and Analytical Determinations of the Stress Systems in a 
Welded Pressure Vessel,” anonymous, published by British Welding Re- 
search Association. Hydrostatic tests of a boiler drum with various 
sizes of nozzle outlets. Stresses around the nozzle openings are mea- 
sured by electrical resistance type strain gages. Strains were measured 
on both the outside and inside surfaces. The tests described in this report 
were carried out during 1946 and 1947 by Babcock and Wilcox, Ltd. 


9. “Analysis of Experimental Data Regarding Certain Design Features of 
Pressure Vessels” by G. J. Schoessow and E. A. Brooks, Transactions of 
the ASME, July, 1950. Analysis of data obtained by strain gage experi- 

mental stress analysis on drum heads and nozzle openings in heads and 

shells. 


“Fatigue Tests of Piping Components” by A. R. C. Markl, Transactions 
of the ASME, Vol. 74, No. 3. Cyclic bending tests on various pipe com- 
ponents such as reducers, curved pipe, miter bends, forged tees and fab- 
ricated branch connections. Results are compared with the strength of 
straight pipe and generalized formulas for the stress intensification fac- 
tors of the piping components are presented. 


11. “Design of Hot Tap Connections in High Pressure Pipelines” by A. J. Del 
Buono, P. L. Vissat and Frank G. Williams, presented at the ASME 
Petroleum Conference, Houston, Texas, September 28-30, 1953. Presents 

a review of the forces acting on branch connections in a pipeline and dis- 

cusses the design of such branch connections. Appendix I summarizes 

some previous research work on branch connections in pressure vessels. 
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Appendix II is a report on tests of an encirclement tee reinforced branch 
connection; hydrostatic tests were run on this reinforced branch connec- 
tion with stresses determined by means of electrical resistance strain 
gages placed on the outside surface. After completion of the strain gage 
tests, the pressure was raised until fracture of the specimen occurred. 


12. “Deformations and Stresses in Circular Cylindrical Shells Caused by Pipe 
Attachments” by N. J. Hoff, J. Kempner, S. V. Nardo, and F. V. Pohle, 
November, 1953, Knolls Atomic Power Laboratory, Report KAPL 921, 922, 
923, 924, 925 and 1025. Theoretical investigation of the stresses in a 
pipeline branch connection under bending loads applied to the attached 
pipe. The load transmitted by the pipe is assumed to be a load acting ona 
line segment slightly shorter than the diameter of the branch pipe. A cor- 
rection for the effects of the finite pipe diameter is obtained in an approxi- 
mate manner by means of the theory of circular plates. 


13. “Researches on Welded Pressure Vessels and Pipelines” by Nicol Gross, 
British Welding Journal, April, 1954. Hydrostatic tests in which strains 
were measured on various types of branch connections, both on the out- 
side and inside surfaces. Covers cyclic pressure tests on various types 
of flat head pressure vessel closures and a forged 6" tee. 


14. “Stress Analysis of Split Tee Welded Sleeve Reinforcement of Welding 
Branch Connection for Hot Tapping,” anonymous, Report by The Ladish 
Company, Metallurgical Department, November 1, 1954. Hydrostatic tests 
of a 24” x 12” branch connection reinforced with a complete encirclement 
saddle. Stresses were determined by means of electrical resistance 
strain gages on both the outside and inside surfaces. Upon completion of 
the strain gage tests, pressure was increased until fracture of the assem- 
bly occurred. 


15. “Stresses from Radial Loads in Cylindrical Pressure Vessels” by P. P. 
Bijlaard, published in Welding Research Supplement of the Welding Jour- 
nal, December, 1954 and December, 1955. Theoretical development of the 
stresses in a cylindrical shell resulting from a distributed radial load of 
rectangular shape and by non-uniform distributed loads giving either cir- 
cumferential or longitudinal localized bending moments on the shell. 
Numerical calculations are given in graphical form and the comparison is 
made with the test data of Schoessow and Kooistra. 


“Welded Tee Connections” by A. G. Barkow and R. A. Huseby, Welding 
Research Council Bulletin Series, No. 22, May, 1955. Hydrostatic tests 
on eight 30” O.D. run x 16" O.D. branch connections. Six of the test 
assemblies were reinforced with saddles, one with a complete encircle- 
ment sleeve, and one with a “Weld-O-let.” Stresses were measured on 
the outside surface by means of electrical resistance strain gages. Upon 
completion of the strain gage test, pressure was increased until fracture 
occurred. 


“Investigation of Static and Fatigue Resistance of Model Pressure 
Vessels” by Julien Dubuc and Georges Welter, The Welding Journal, 

July, 1956. Hydrostatic and cyclic pressure tests of 13 1/2 " O.D. x 3/4 
wall pressure vessels. Vessels had 2" O.D. x 3/8 ” wall nozzles; initial 
failures in cyclic pressure tests occurred at these nozzles. Notches and 
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grooves, machined in the outside surface of the shell to simulate acciden- 
tial notches that might be present in an actual pressure vessel, proved to 
be lower stress raisers than the nozzle outlets. Two materials were in- 
vestigated, one with about 35,000 yield strength, 56,000 ultimate strength; 
the other with about 77,000 yield strength, 96,000 ultimate strength. 
Cyclic pressures ranged from zero to approximately the yield pressure 
of the vessels. 


con 
j 
| 
4 
4 
4 
a 
q 
Migs. 
| 
J 
3 
q 
at 
- 
4 
4 
a 
q 


Paper 1194 


Journal of the 
PIPELINE DIVISION 


Proceedings of the American Society of Civil Engineers 


FLOW OF NATURAL GAS IN PIPELINES* 


William T. Ivey,! M. ASCE, and James H. Dorough2 
(Proc. Paper 1194) 


SYNOPSIS 


Consideration is given to the various flow formulae in use in the design of 
natural gas transmission pipelines, and to the performance ratings of these 
lines. Generally, all formulae are derived from the general energy equation, 
the basic difference being found in the methods of considering the friction co- 
efficient, and in the adjustment for the deviation of gases for supercompres- 
sibility. Each formula, after applying all adjustments, assumes the incor- 
poration of an “efficiency factor” to record the deviation of the formula from 
a selected norm, although no basic norm, or perfect operating condition, has 
been determined by the industry. 

The need for such a norm in the form of a standard formula based on maxi- 
mum possible efficiency of operation, against which actual line performance 


The merits of a proposed form for this standard equation is presented 
herein. Each of the basic differences is examined in order to attempt to 
simplify the mechanical application of the formula, to permit its more rapid 
use by manual application and by digital computers, thus permitting a more 
rapid and thorough examination of the many problems presented by the in- 
creasingly complicated pipeline networks of the present day. 


HISTORICAL 


The history of the development of present day formulae for flow of gas in 
pipelines had its beginning around the turn of the century. The earlier 


Note: Discussion open until August 1, 1957. Paper 1194 is part of the copyrighted 
Journal of the Pipelines Division of the American Society of Civil Engineers, Vol. 
83, No. PL 1, March, 1957. 


a. Presented at a meeting of the Pipeline Division, ASCE, February 18, 1957, 
Jackson, Miss. 

1. Chf. Statistician and Planning Engr., Southern Natural Gas Co., 
Birmingham, Ala. 

2. Planning Engr., Southern Natural Gas Co., Birmingham, Ala. 
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applications did not consider deviation, and generally applied a constant fric- 
tion coefficient. Such formulae were the formulas of Cox, proposed in 1902, 
and Rix, proposed in 1904. Other early formulae expressed the friction co- 
efficient as a function of the diameter. Such formulae were those of Oliphant 
(1904), Unwin (1904), Weymouth (1912), California (1925) and Spitzglass 
(1926). A third group of formulae expressed the friction coefficient as a 
function of Reynolds’ criterion, a factor considering diameter, velocity, den- 
sity, and viscosity. Such formulae were: Fritzsche (1908), Lee (1914), White 
(1929) and McAdams and Sherwood (1926). The Panhandle formulae fall in 
this category. The Panhandle formulae make the friction coefficient a func- 
tion of Reynolds’ number, with a constant viscosity factor. 

The first attempt to produce a standard formula was made in 1935 when 
T. W. Johnson and W. B. Berwald presented a joint report on flow of natural 
gas through high-pressure transmission lines, under the sponsorship of the 
Bureau of Mines and the Natural Gas Department of the American Gas Associ- 
ation, issued as Monograph No. 6 of the Bureau of Mines. Numerous data on 
gas flow through commercial transmission lines were obtained and analyzed. 
This was probably the most thorough and comprehensive analysis of flow of 
gas in pipelines to that date. Johnson and Berwald concluded that for the 
larger diameter lines, free from condensates and other foreign materials, 
and operating under steady flow conditions, the metered rates of delivery 
agreed more closely with the rates calculated from Weymouth’s formula than 
with those calculated from any of the other formulae; and in almost every 
test, under these conditions, the volume calculated from Weymouth’s formula 
was within a few percent of the metered delivery. Further experiments, un- 
der actual conditions of natural gas transmission, were recommended. 

The work of Johnson and Berwald was extended in 1956 by R. V. Smith, 
J. S. Miller, and J. W. Ferguson, in a joint report on flow of natural gas 
through experimental pipelines and transmission lines. Issued under the 
sponsorship of the Bureau of Mines and the Pipeline Research Committee of 
the American Gas Association this publication appeared as Monograph No. 9. 
The report was intended primarily to present the results of tests conducted 
on experimental pipelines, and to apply these data to a further improvement 
of the formula. Tests were made on experimental pipelines constructed of 
2, 4, 6 and 8 inch seamless pipe, 8 inch welded pipe, and a specially prepared 
3 inch superfinished pipe. 

The friction coefficient, “f” was redefined as “resistance coefficient,” and 
a new factor “V1/f” termed “transmission factor” was introduced. 

Transmission factors for the various experimental pipelines were deter- 
mined for 2, 4, 6 and 8 inch pipe. Relationships were established between 
transmission factors and the pertinent variables: Reynolds number, absolute 
roughness, and internal radius of the pipe. Reynolds numbers for these de- 
terminations of transmission factors ranged from a low of 49,900 to 
36,000,000. At Reynolds’ numbers ranging from 250,000 to approximately 
850,000 transmission factors for the two inch seamless were found to be re- 
lated to Reynolds’ numbers, internal radius of the pipe, and the absolute 
roughness of the internal surface of the pipe. At Reynolds numbers above 
850,000 the transmission factors were constant and independent of the 
Reynolds’ number, being influenced only by the absolute roughness and the in- 
ternal radius of the pipe. For larger diameters, extrapolated data indicated 
the break between transition and fully rough flow moved to higher Reynolds’ 
numbers. 
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Analysis of pipeline flow efficiency field test data by the Institute of Gas 
Technology in 1952 and 1953 disclosed inconsistencies in the data resulting 
from differences in measuring technique. To satisfy the need for proper 
evaluation of field measurement procedures, and for development of criteria 
for judging the validity of field data, a set of preliminary recommendations 
for efficiency testing was developed in 1954. In May, 1956, a formal pro- 
cedure for pipeline efficiency testing was presented by R. F. Bukacek and 
R. T. Ellington. 

Currently the American Gas Association, working with the Institute of 
Gas Technology at Chicago, through its Project NB-13, “Pipeline Efficiency 
Flow Standards” is in the process of analyzing an appreciable quantity of 
test data. These data, provided through the cooperation of a number of pipe- 
line companies, cover field tests of pipeline installations and include a wide 
range of diameters and various flow rates. No conclusions have, as yet, 
been published by this committee. 

It is not proposed here to anticipate the conclusions of this committee. It 
is desired to consider certain data, consisting of over two hundred tests 
made on pipelines of Southern Natural Gas Company, and to present certain 
conclusions reached upon examination of that data. The tests were made un- 
der the supervision of the Operating Department of Southern Natural Gas 
Company. 


The Basic Formula 


Much information is available in the literature on flow of gas in pipelines. 
Numerous varying formulae are offered for the practical solution of the 


problem. They all originate from the basic flow equation, as derived from 
the general energy equation. The basic derivation of the general flow formu- 
la appears in the work of Johnson and Berwald in U. S. Bureau of Mines 
Monograph No. 6 and will not be repeated here. 


The general flow formula: 
2 2 1/2 
- Pp) a 
GTLf 


= volume of gas in cubic feet per hour at pressure base Po and tem- 
perature base To; 

= numerical constant 1.6156; 

= temperature base, defining one cubic foot of gas, degrees F absolute; 

= pressure base defining one cubic foot of gas, pounds per square inch 
absolute; 

= inlet pressure, pounds per square inch, absolute; 

= outlet pressure, pounds per square inch, absolute; 

= internal diameter of pipe, inches; 

= specific gravity of gas (air = 1.000); 

= temperature of flowing gas, degrees F absolute; 

= length of pipe in miles; 

= coefficient of friction. 
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This gives no recognition to the deviation factor, which is recognized by 
Smith, Miller, and Ferguson in Monograph No. 9 as follows: 


° GTL 


where 


f = resistance coefficient, dimensionless; 
V 1/f = transmission factor, dimensionless; 
Z = average effective compressibility factor of gas, dimensionless; 


It is noted that “f,” formerly defined as “coefficient of friction,” is now 
defined “resistance coefficient,” and a new term “ V 1/f ,” defined as “trans- 
mission factor” is introduced. This change was made in an effort to estab- 
lish clear and concise terminology for these terms. 

This formula is the basic flow equation from which all other formulae 
are derived. At this point, theory gives way to laboratory or field test in the 
evaluation of the “f” term, or the term that will bring the equation in balance. 
Three basic approaches have been used in the determination of the various 
formulae proposed: (1) the application of the friction factor, or resistance 
coefficient as a constant; (2) its application as a function of the diameter of 
the pipe only; and (3) as a function of the Reynolds number. 


Reynolds Number 


The Reynolds’ criterion, a dimensionless ratio commonly used to charac- 
terize the conditions of fluid flow in circular pipes, is stated as: 


_ DUS 


where: 


R = Reynolds’ criterion 
D = Internal diameter of pipe 
U = Average linear velocity 
S = Density of fluid 

V = Viscosity 


Reynolds observed that at low velocities flow proceeded in a steady 

i straight-line manner, in which it appeared that all of the particles moved 

; Parallel to the walls of the tube. As velocities increased, a certain critical 
change in the character of the flow occurred, and turbulence appeared. 
Smith, Miller and Ferguson noted a change in the relationship between the 
transmission factor and Reynolds’ number as this critical point was reached. 
They stated that transmission factor variations, for the experimental data on 
two inch diameter pipe, depended on Reynolds’ numbers only for ranges in 
Reynolds’ numbers from 49,900 to 250,000. At Reynolds numbers ranging 
from 250,000 to approximately 850,000 transmission factors were related to 
Reynolds’ numbers, internal radius of the pipe, and the absolute roughness of 
the internal surface of the pipe. At Reynolds’ numbers above 850,000, for the 
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pipes considered, transmission factors were constant and independent of 
Reynolds’ numbers. They were influenced only by the absolute roughness and 
the internal diameter of the pipe. For larger diameters, extrapolated data 
indicated the break between transition and fully rough flow moved to higher 
Reynolds’ numbers. Smith, Miller and Ferguson interpreted these three 
areas of flow as “smooth,” “transition,” and “fully rough,” which compare 
with Reynolds observations of “steady” and “turbulent” areas of flow. 


Consideration of Test Data 


Over the past several years, a continuing series of careful tests have been 
carried out on the pipelines of Southern Natural Gas Company under the 
supervision of the Operating Department of Southern Natural Gas Company. 
Pipelines tested consisted of sizes ranging from ten-inch to twenty-four inch, 
inclusive. Flow measurement was made by metering where possible, and 
checked by means of ammonia tracer. Due to the variable distances of ori- 
fice measurement from the point of test, the ammonia tracer method was 
found to give more consistent results. These tests are summarized in 
Appendix A. 

Volumes of gas are stated in Mcf per day at a pressure base of 14.73 psia 
and are as measured by ammonia tracer. Where ammonia tracer measure- 
ment was not available, metered flow is used and the values so obtained are 
set out by an asterisk. The table shows the measured flow, the internal di- 
ameter of the pipeline in inches, the difference between the squares of the 
initial and terminal absolute pressures, the average absolute pressure in the 
test section, the gravity of the gas, the absolute temperature, the deviation 
factor and the gas viscosity. An adjustment is made for elevation. This ad- 
justment is discussed in Appendix B. 

A solution is made for Reynolds’ number, resistance coefficient, and 
transmission factor. 

Certain conclusions are immediately evident: 


 e Two hundred and thirty-three tests were conducted in actual operation 
of a large transmission system. Under varying ranges of normal operat- 
ing conditions, results were obtained from varying lengths of pipe ranging 
in nominal size from ten to twenty-four inches and under pressures rang- 
ing from three hundred to twelve hundred pounds per square inch gage. 
During this entire series Reynolds’ numbers varied from a low of 
2,202,511 for one twelve inch pipe to a high of 13,272,727 for one eighteen 
inch pipe. That is, all of the tests made showed Reynolds’ numbers in 
ranges much higher than Reynolds’ numbers in the “smooth” and “transi- 
tion” ranges defined by Smith, Miller, and Ferguson. They were also well 
above the Reynolds’ number which, according to Smith, Miller and 
Ferguson, define the limit above which the transmission factors are in- 
fluenced only by the absolute roughness and internal diameter. 

Therefore, if this conclusion can be borne out by further tests, normal 

pipeline operation need not be concerned with the Reynolds’ criterion, but 
only by diameter and by absolute roughness. 


2. Plots made of transmission factor against Reynolds’ number for 16, 18 
and 24 inch diameters from the test data available show a broad band of 
points at each Reynolds number and no clearly defined relationship between 
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Reynolds’ number and transmission factor within the range of the tests 
made. 

If the influence of Reynolds’ number is not felt in rough, turbulent flow, 
but rather if the transmission factor is influenced only by the absolute 
roughness and internal diameter, then the data should be examined by 
sizes of pipe, and roughness factors calculated for each pipe size. The 
formula used by Smith, Miller and Ferguson is their resistance equation 
(Formula 28, Monograph 9): 


V 1/f = 4 log. iar (4) 


where r = internal radius of pipe in inches 
k = absolute roughness coefficient 


The results of this examination are tabulated in Table I. 

Each line in Table I represents a group of tests on a single line or 
lines of the same diameter and in the same continuing run of pipeline. 

The nominal pipe size is shown in the first column as a matter of conveni- 
ence, but all calculations are based on the internal diameter shown in the 
second column. The number of tests made, and the maximum, minimum 
and average transmission factor are shown for each group of tests. The 
absolute roughness factor is calculated by the equation (4), on the assump- 
tion that all Reynolds numbers exceed the critical criteria for rough tur- 
bulent flow, and therefore, the transmission factor is a function of diame- 
ter and absolute roughness only. 

Line 1, Table I, gives the results of a total of twenty tests on a new 
solid welded twenty-four inch pipeline operating in the range of 500 to 750 
pounds per square inch gage. The transmission factor varies from 19.741 
to 20.933 with an average transmission factor of 20.195, or for a deviation 
of less than 3%, plus or minus. The calculated absolute roughness coeffi- 
cient is 0.000777. This line is one of the most recent installations and has 
been in operation for the entire period on clean gas. Smith, Miller and 
Ferguson reported absolute roughness from 0.000553 to 0.00185, and used 
an absolute roughness factor of 0.0007 in the calculation of the limits of 
the transition region. 

For the purpose of this analysis, it may be assumed that an absolute 
roughness of 0.000777 and a transmission factor of 20.000 for 24 inch pipe 
indicates a line operating at high efficiency. 

Line 2 gives the results of eleven tests on twenty-two inch pipe and 
shows an absolute roughness factor of 0.000721. Line 5 gives the results 
of 18 tests on twenty-two inch pipe with an absolute roughness of 0.000624. 
Line 17 gives the results of four tests on eighteen inch pipe with an abso- 
lute roughness of 0.000691. Line 23 gives five tests on sixteen inch with 
an absolute roughness of 0.000792 and line 24 gives five tests on eighteen 
inch with 0.000789. Line 33 gives 17 tests on 18 inch with 0.000761. All 
of these tests may be interpreted as tests on pipe in good condition, 
operating efficiently. 

Line 3 gives the results of eight tests on twenty-inch pipeline with an 
absolute roughness of 0.002596. Line 4 gives the results of six tests on 
twenty-four inch with an absolute roughness of 0.001503. Line 7 gives four 
tests on twelve-inch with 0.000848 coefficient, line 8 gives sixteen tests on 
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TABLE I 
COMPARISON OF TRANSMISSION FACTORS AND ABSOLUTE 
ROUGHNESS COEFFICIENTS BY PIPE SIZES FROM TESTS 
MADE BY SOUTHERN NATURAL GAS COMPANY 


PIPE SIZE Absolute 
INCHES Number TRANSMISSION FACTORS Roughness 
Line 0O.D. ee of Tests Minimum Maximum Average Coefficient 


23.500 19.741 20.195 0.000777 
21.375 19.630 20.160 0.000721 
19.313 16.265 17.591 0.002596 
23.188 17.565 19.026 0.001503 
21.375 19.245 20.413 0.000624 
17.375 18.315 19.350 0.000935 
12.125 18.257 18.893 0.000848 
19.313 17.678 18.921 0.001329 
23.188 16.860 19.207 0.001354 
13.375 15.811 0.005518 
15.375 18.257 0.001552 
12.125 16.514 0.003372 
19.313 16.566 18.846 0.001389 
19.313 18.474 19.351 0.001038 
13.375 17.413 18.752 
15.438 18.788 19.400 
17-375 19.830 19.875 
15.438 19.2h5 
17.375 19.181 19.355 
15.375 18.380 19.211 
12.250 18.898 19.503 
10.250 18.896 
15.438 19.432 
17.375 19.645 
19.250 17.230 
21.375 19.288 
17.500 18.178 
21.375 18.709 
17.375 18.768 
21.375 19.415 
23.390 20.000 
23.500 18.688 0.001850 
17.375 19.706 0.000761 
19.313 17.639 18.620 0.001581 
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twenty-inch with 0.001329, and line 9 gives six tests on twenty-four inch 
with 0.001354. All of these lines are known to have contained foreign mat- 
ter, presumably liquids, at the time of the tests. 

Lines 10, 11 and 12 were tests on field gathering lines which from time 
to time had received liquids. These tests may be presumed to show re- 
sults on lines which have operated at poor efficiency. 

Lines 25, 26, 27, 28, 29 and 30 record tests on the oldest portion of the 
system, a great part of which consisted of coupled lines. These tests may 
be presumed to show results on old lines which have operated for a num- 
ber of years. 

Thus, a range of absolute roughness factors may be determined to show 
operating conditions from maximum efficiency to poor efficiency. 

It is concluded that, for the tests considered, the condition of the pipe- 
line, with respect to liquids or other foreign matter, and with respect to 
the age of the pipeline and type of construction, inclusion of couplings, 
valves, and cross connections, has a much greater influence on the trans- 
mission factor than the pipe diameter or Reynolds’ criteria. 


The Dilemma of the Gas Industry 


It is apparent from the above discussion that it is not possible to specifi- 
cally express the percent efficiency of a pipeline without a basic norm or 
standard base of reference. Each gas transmission company rates its pipe- 
lines according to its own formula, assuming friction coefficients based on 
either Reynolds’ numbers or pipe diameters, and on roughness coefficients 
related or not related to diameters. A line that may be rated at high effi- 
ciency by one formula will be rated at a low efficiency by another. In fact, by 
some formulae, a line may be rated at over 100% efficiency. 

It is evident that line efficiency cannot be adequately determined by fric- 
tion coefficients based on a constant without reference to line size. It also 
appears that in many actual operating conditions the friction coefficient or 
transmission factor is independent of Reynolds’ number. May it be assumed 
that a practical solution for commercial pipelines under normal operating 
conditions is in the area of high Reynolds’ numbers with transmission factors 
based on diameters and absolute roughness alone? 

Smith, Miller and Ferguson suggest a tentative average absolute roughness 
of 0.0007, and agree that for fully rough flow the transmission factor is not a 
function of Reynolds number. On this basis, transmission factors are calcu- 
lated for varying pipe diameters in Table I. 

Such transmission factors are not offered as final, but are suggested as 
probably representative of transmission factors for turbulent flow in areas 
of high Reynolds’ numbers. Acceptance of a basic absolute roughness coef- 
ficient for pipe of maximum practical smoothness will permit the determina- 
tion of transmission factors which might be used as standards of maximum 
line efficiency. 


The Deviation Factor 


It has been shown that natural gas does not follow the ideal gas relation- 
ship, but that, at high pressures greater volumes of gas may be compressed 
into a given space than the ideal gas laws indicate. Failure to recognize this 
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TABLE II 
TRANSMISSION FACTORS FOR ROUGH TURBULENT FLOW 
THROUGH PIPELINES OF VARIOUS DIAMETERS 


| Vi7F = k log d+ 14.9; K = 0.0007 ] 


Pipe Size; Inches Transmission 
0.D. Wall I.D. Factor 


4 1/2 17.362 
6 5/8 18.048 
8 5/8 18.539 
10 3/4 10.250 18.943 
12 3/4 12.250 19.253 
14 13.500 19.421 
16 15.375 19.647 
18 17-375 19.860 
20 19.3125 20.043 
22 21.250 20.209 
ok 23.1875 20.361 
26 25.125 20.500 
30 29.125 20.757 


x 
x 
x 
x 
x 
x 
x 
x 
x 
x 
x 
x 
x 


deviation in the past has not introduced a material error because of the low 
pressures at which earlier pipelines operated. However, with continuing in- 
creases in operating pressures, recognition of this deviation or supercom- 
pressibility characteristic of the gas becomes more and more important. 

This subject has been most thoroughly covered by Dr. George Granger 
Brown in his paper, “Deviation of Natural Gas from Ideal Gas Laws.” Dr. 
Brown points out that although numerous equations of state have been pro- 
posed, the most convenient method for computing the pressure-volume- 
temperature relationship of natural gas is by the use of the so-called com- 
pressibility factor, which represents the deviation of the gas in question from 
the ideal gas laws, as expressed by the equation 


PV = ZNRT (5) 


where 


P = pressure in pounds per square inch absolute 
V = Volume in cubic feet 

Z = the compressibility factor 

N = number of pounds mol 

R = 10.71 for all gases 

T = absolute temperature 


The compressibility factor Z is a dimensionless intensive factor 
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independent of the extent or weight of the gas, and determined by the charac- 
ter of the gas, the temperature, and the pressure. Once Z is known or deter- 
mined, the calculation of the pressure-temperature-volume relationship may 
be made with as much ease at high pressure as at low pressure. 

According to the theorem of corresponding states, the deviation of an ac- 
tual gas from the ideal gas law is the same for different gases when at the 
Same corresponding state. The same corresponding states are found at the 
same fraction of the absolute critical temperature and pressure which are 
known respectively as 


Reduced Temperature . = (6) 


(7) 


Reduced pressure 


where 


Tc = the absolute critical temperature 
P,. = the absolute critical pressure 

T =the absolute temperature at which the gas exists 
P =the absolute pressure at which the gas exists 
Ty = reduced temperature 

P, = reduced pressure 


Critical temperature is defined to be the temperature above which a pure 
gas cannot be liquefied by pressure alone. Critical pressure is defined to be 
the pressure below which a pure gas may exist in a gaseous state in equilib- 
rium with its liquid. 

Therefore, if the theorem of corresponding states can be applied without 
appreciable error, all gases would have the same value of Z at the same re- 
duced temperature and pressure, and a plot of Z for methane, as a function 
of reduced temperature and pressure, can be applied to determine the un- 
known value of Z for some other gas if we know or can determine the critical 
temperature and pressure of the second gas. 

In the previously mentioned paper Dr. Brown presented such a plot for 
methane and for natural gas which has found wide use in the industry in the 
determination of compressibility factors. 

Based on this chart, deviation factors are determined for each of the vari- 
ous gas mixtures found in the system on which gas tests were made. Charac- 
teristics of the gas considered are determined by laboratory analysis. Indi- 
vidual absolute critical temperatures and pressures of each constituent are 
set down and the pseudo-reduced temperature for the mol-fraction calculated. 
The arithmetic total of these fractions will give the absolute pseudo-critical 
temperature and absolute pseudo-critical pressure for the gas mixture con- 
sidered. From these values the pseudo-reduced temperatures and pseudo- 
reduced pressures for varying actual temperatures and pressures are deter- 
mined. From these factors, the deviation factors may be read from the chart. 

As Dr. Brown points out, like all correction factors, the value of Z 
depends upon experimental data for its determination. 

The standard apparatus for experimentally determining supercompressi- 
bility is either the Bean or the Burnette. The Bean apparatus consists of a 
high pressure gas bomb of known volume and a calibrated burette for measur- 
ing the volume of the gas at atmospheric pressure. It is simply an adaptation 
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of the classic laboratory procedure of expanding an unknown quantity of gas 
in successive increments from a high pressure bomb into a calibrated 
burette at essentially atmospheric pressure, and computing the total volume 
from a summation of these increments. The operation is tedious; mistakes 
are easy to make, and one mistake can ruin the whole run. The Burnette ap- 
paratus consists of two adjacent high pressure chambers of a known volume 
ratio, the gas being expanded from the first chamber into the evacuated 
second chamber, with successive expansions being made after re-evacuation 
of the second chamber. The deviation of the gas is computed from a compari- 
son of the volume ratio of the chambers to the pressure ratio before and af- 
ter the expansion. The operation is very easy and not conducive to mistakes. 
The methods are of comparative accuracy, agreeing to within about 0.1%. 
However, due to ease of operation, the Burnette is most used. 

In 1955, the American Gas Association, through PAR Project NX-7, pre- 
pared such data based on the specific gravity method of super-compressibili- 
ty factor determination. Tables of supercompressibility factors Fpy for 
natural gas were issued for nominal temperatures from 0° to 180° F, for 
specific gravity from 0.554 to 0.750, and for pressures of from 0 to 3,000 
psig. Separate tables were included for the determination of supercompres- 
sibility factors for natural gas containing nitrogen and or carbon dioxide. 


By definition F BRT (8) 


From Equation 5 


(9) 


Therefore, 
(10) 


In Table I, a comparison is made between Z factors calculated by the 
Granger Brown method and Fpy factors as taken from the A.G.A. tables. 
The differences are negligible, and it may be concluded that deviation factors 
may be determined by either method. 


A Proposed Formula for Flow of Natural Gas 


It appears, at least for the data available here, that the following con- 
clusions can be drawn: 


For turbulent flow the transmission factor is not a function of Reynolds 
number 

Normal commercial operation will fall within this range 

The transmission factor will be a function of diameters and absolute 
roughness 

The absolute roughness coefficient for pipelines in good condition 
should approximate 0.0007 

The deviation factor may be obtained by the George Granger Brown 
method or from the A.G.A. tables of deviation 

All of these factors may readily be incorporated in the basic flow 
equation. 
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(1) 


TABLE III 
COMPARISON OF DEVIATION FACTOR Z OF WESTFIELD GAS 


PL1 


CALCULATED BY GRANGER BROWN METHOD 


WITH F 


(2) 


Pressure wen 
Psig at 60°F 
100 0.983 
200 0.967 
300 0.950 
400 0.934 
500 0.918 
600 0.903 
700 0.890 
800 0.876 
900 0.862 
1000 0.8h8 
1100 0.837 
1200 0.825 


Proposed Formula 


Q=k Vi/t V1/z 


FR 


L 


(3) 
1 


7, 
at 60°F 


1.0173 
1.0341 
1.0526 
1.0707 
1.0893 
1.1074 
1.1236 
1.1416 
1.1601 
1.1792 
1.1947 
1.2120 


OM A.G.A. TABLES 


(4) 


Vin 


1.0086 
1.0169 
1.0260 
1.0347 
1.0437 
1.0523 
1.0600 
1.0685 


GTL 


1.0859 
1.0930 
1.1009 


Therefore, for a design formula applicable to normal commercial use of 
pipelines in good condition, the basic flow equation may be used as follows: 


2 2. .5 i/2 
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(5) 


A.G.A 
Tables 
_ 
1.0083 
1.0164 
1.0246 
1.0329 
1.0416 
1.0501 
1.0589 
1.0677 
1.0763 
1.0847 


1.0933 
1.1015 


(11) 


For normal use, a constant (Kj) may be determined to include all constant 
factors, including diameter and transmission factor and a table prepared to 
give K; for each diameter. The equation then becomes: 


1/2 


2 
Q= 


(12) 


At this point consideration is given to the application of the factor V 1/Z. 
In most applications, this is taken to be the deviation for the gas at the mean 
effective pressure between the initial and terminal pressure. 
A compressibility factor may be determined, as more fully outlined in 
Appendix B, which may be applied directly to the pressure. 
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The modified flow formula then becomes: 


E - (Po Ve (13) 


L 


A plot of deviation factors for a representative gas may then be deter- 
mined for varying pressures, and a tabulation prepared adjusting the actual 
pressure to an adjusted pressure recognizing the deviation factor. 

In this form, through the use of a table of constants taking into account 
diameter and transmission factor for each pipe size, another a table giving 
squares of absolute pressures, adjusted for the deviation of a typical gas for 
varying gage pressures, the formula is readily applied without the use of 
graphs, charts, logarithms, and without cut andtry. It may be handled rapid- 
ly in manual operation, with the use of a desk calculator, or may be keyed in- 
to the operation of any of the basic digital computers. 

In this manner, through the use of digital computers, it will be possible to 
make a much more rapid and more thorough examination of the many and 
complicated problems both of construction and of operation facing the in- 
creasingly complicated pipeline networks of today. 


CONCLUSIONS 


Many different formulae have been presented for determination of flow of 
gas in pipelines. All are derived from the basic flow equation. Chief differ- 
ence is in the determination of friction coefficient. Some formulae set a 
basic friction factor for all conditions; some relate friction factor or trans- 
mission factor to diameter and absolute roughness. Others relate transmis- 
sion factor to Reynolds number. Early formulae did not adjust for deviation. 
It is generally accepted today that such adjustment is necessary. 

For practical application, based on the tests submitted herein, it is rea- 
sonable to conclude transmission factor is not a function of Reynolds’ number 
for turbulent flow in the area of normal commercial operation, but that trans- 
mission factor is a function of diameter and absolute roughness. It appears 
that variations of absolute roughness due to varying pipe conditions is of 
greater significance than pipe diameter. It is urged that a basic absolute 
roughness factor, such as 0.0007, be agreed upon as being the norm against 
which pipeline efficiency can be determined. 
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APPENDIX B 


Correction for Elevation Change 


In areas in which there are material differences in elevation it is neces- 
sary to adjust for these differences. Such adjustments have been made in the 
analysis of test data given in Appendix A. This correction has been made by 
the application of an adjusting factor, shown here as an adjustment to the dif- 
ferences of the squares of the absolute pressure in Equation 12. 


(P,2 0.0375 Hy) _Pavg x\/1/ (14) 
L 


Q= Kk) 


where the elevation factor term is 


0.0375 G (H, ~ H,) Pavg” 
Zavg . Tavg 
G = gravity of flowing gas (air = 1.000) 
Z = average supercompressibility factor 
T = average flowing temperature °F absolute 
P = average flowing pressure psia 


H, = elevation at inlet of pipe section (feet) 
Ho = elevation at outlet of pipe section (feet) 


Calculation for Reynolds’ Number 


-4 9G 
R = 135.321 x10 x Vd 


where 


Q = Mcfd at 60 °F and 14.73 psia 

G = gravity of gas (air = 1.000) 
*V = viscosity of gas (lb. per sec. ft.) 
d = internal diameter of pipe inches 


Application of Deviation Factor 


Beginning with the proposed formula (Equation 12) 


Q=K x 1/z (12) 
2 L 


This formula expresses the deviation as a function of the mean effective 
pressure between the initial and terminal pressure of the span being: 


* V in Appendix A shown times 108, 
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considered. This formula may also be expressed 
- m2) a2] (16) 
Q = Ky 


The literature suggests that 


March, 1957 


1/z = (1 + JPy) (17) 


where JPy is the average deviation factor at the mean effective pressure 
and may be expressed as (1/Z - 1) for the pressure considered. Py, the 
mean effective pressure, can be equated to P1 and P92 as follows. 


|= (18) 
Foe 


Substituting Equation 17 in Equation 16 gives 


1/2 
(1) 


Substituting Equation 18 for Py in Equation 19 reduces to 


1¢ 


- Pp3 1/2 
Q=K, (Py2 - Po2) + J 2/3 ( ) - 


L 


Simplifying and grouping 


_— 1/2 


L 
- + J 2/3 - 2/3 1/2 (22) 
L 
1/2 
Q=K, (1+5 2/3 Pl) - (1+ 2/3 P 
(23) 


Then, expressing 


(1 + J 2/3 Py) = »,° (24) 


Where b is a factor for the determination of the deviation of the gas at the 


: 
: 
q 
A 
20) 
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particular pressure being considered, the formula becomes: 


L 


which formula is formula 13 above mentioned. 

Thus, a table may be prepared making the appropriate adjustment for 
deviation of the gas considered, and relating the square of the absolute pres- 
sure adjusted for deviation to the corresponding gage pressure. This table 
will give factors which can be substituted in the proposed equation which will 
adjust for deviation. 


APPENDIX C 
Verification of Equation 15 


“The Reynolds number, a dimensionless ratio commonly used to charac- 
terize the conditions of fluid flow in circular pipes, is given in the literature 
as: 


p 


(26) 


ae 


where 


R = Reynolds number; 

D = Diameter of pipe; 

u = average linear velocity; 

p = density of fluid; 

uM = dynamic viscosity of fluid; 

v = kinematic viscosity of fluid (U/p). 


The above relationship has been expressed by Biddison! in natural-gas 
engineering terms as 


R = 0.011459 SCP, (27) 


vaT,, 


where 


Q = rate of flow, cubic feet per hour at a pressure and temperature base 
of Ph p.s.i.a. and Tp, OF abs.; 

G = specific gravity of gas (sp. gr. of air = 1.0); 

V = viscosity of gas, pounds per second foot; 

d = internal diameter of pipe, inches 


The constant in equation (26) adjusts for the mixed units used in the 
equation. ”2 


1. Biddison, P. M. Gas Flow Computations. Proc. Am. Gas Assoc. Nat. Gas 
Dept. 1941, pp. 51-88. 
2. Bureau of Mines Monograph Number 9, Page 32 (Quotation taken from) 


rity 


- 
a 
= 
ae 
pa 
a 
a 
Son 
| 
| 
| 
om 
a 
al 
= 


1194-28 PL1 March, 1957 


Equation 15 in our paper is obtained by substituting in Equation 26 above 
the following 


1000 


P, = 14.73 Psia 
T, = (60° F + 460) = 520° abs. 


where 


Q, = rate of flow, thousand cubic feet per day (Mcfd) at a pressure and 
temperature base of 14.73 psia and 60° F abs. 


Equation 26 becomes 


14.73 1000 
R = 0.011459 x 590 


Then 


R = 135.321 x 


(Equation 15) 
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Journal of the 
PIPELINE DIVISION 


Proceedings of the American Society of Civil Engineers 


PREDICTION OF SURGE PRESSURES IN OIL PIPELINES* 


A.M. ASCE, and Edwin J. Waller,? J.M. ASCE 


(Proc. Paper 1195) 


Robert D. Kersten, 1 


INTRODUCTION 


Pressure surges have long been the bane of the pipeliner’s existence. In 
spite of the advances in technology and equipment, surge problems for the 
most part have remained unsolved and are generally regarded as a neces- 
Sary evil. undoubtedly, this has been so because the various water hammer 
theories were not fundamental enough in their approach to permit wide- 
spread application to a variety of pumping situations. In addition, many com- 
mercial desurgers or pulsation dampeners while helpful in some cases are 
frequently detrimental in others. 

The net result is that the pipeliner is still concerned with: 


1. The inability to predict the points of extreme fluctuating pressure ina 
pumping system, thereby resulting in the overdesign of the entire 
system. 

2. Breakdowns of the pumping system from elastic failure of the pipes 
and pumps due to extreme fluctuating pressures. 

3. Loss of revenue becayse “surges” necessitate the establishment of a 
line limit below the normal rating of the equipment. 


It is the purpose of this paper to outline a method of analysis based on a 
fundamental mathematical approach and to demonstrate the validity of this 
method through application to an actual field-scale pumping pipeline system. 
The future applications of this method should materially aid the designer in 
predicting points of extreme fluctuating pressure, resulting in better over- 
all designs and appreciable monetary savings. 

Reduced to its most basic concept, pipelining is transportation. Trans- 
portation is movement, and moving objects have inertia and resistance. 


Note: Discussion open until August 1, 1957. Paper 1195 is part of the copyrighted 

Journal of the Pipelines Division of the American Society of Civil Engineers, Vol. 

83, No. PL 1, March, 1957. 

a. Presented at a meeting of the Pipeline Division, ASCE, February 18, 1957, 
Jackson, Miss. 

1. Research Engr., The Carter Oil Co., Tulsa, Okla. 

2. Associate Prof. of Civ. Eng., Oklahoma Inst. of Technology, Oklahoma 

Ag’n. of Mech. College, Stillwater, Okla. 
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Therefore, each investigator in this field has succeeded only to the extent to 
which he has applied fundamental engineering techniques. 


Development of a Method of Analysis 


Previous Investigations 


The earliest noteworthy contributions to the water hammer or pressure 
surge theory were those of Allievi(1) and Joukowski.(2) Joukowski postulated 
the classic water hammer formula for the increase in pressure due to the 
extinguishing of a given velocity of flow. The work of Allievi represented the 
first fundamental mathematical approach to variable flow. Although the basic 
differential equations of variable flow were not explicitly solved for the 
steady-state case until nearly fifty years later, the importance of Allievi’s 
work should not be overlooked. 

Many of the water hammer and pressure surge studies worthy of note 
draw heavily on the work of Allievi and Joukowski. The most significant con- 
tributions in this field have been summarized by the American Society of 
Mechanical Engineers Symposium on Water Hammer(3) and by Wyatt 4) and 
will not be repeated here. It should be noted, however, that for the most part 
the methods of analysis outlined in these papers were of a semi-empirical 
nature and lean heavily on graphical techniques. 


Results of Recent Research 


Mr. Waller(5,6) proposed a solution to the fundamental equations of vari- 
able flow postulated by Allievi,(1) Rayleigh,(7) Lamb,(8) and others. The 
basic differential equations found in this analysis were of the same form as 
those found in the transmission of energy by heat, sound, and electricity. It 
was natural, therefore, that use should be made of the mathematical tool of 
wave analysis. In the process of solution of the differential equations, the 
assumption was made that the principle of superposition was valid for this 
case. The resulting solution consisted of a transient term and a steady-state 
term. 

The transient term was found to be the water hammer solution and obvious- 
ly damps out after a period of time. Analyzing the steady-state term inde- 
pendently, the following result was obtained for the fluctuating pressure (P) 
and the alternating volumetric flow rate or volume current (Q) at any point 
“x” distance from the receiving end: 


Note: Letter symbols adopted for use in this paper are completely 
defined and arranged in alphabetical order for convenience of 
reference in the appendix. 


(1) 


Py, = Py Cosh¥ x + SinhY¥x 


P. 
= Cosh ¥ x + SinhY x (2) 


where Z, and Y are constants defined as: 
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(3) 
(4) 


where a is the attenuation coefficient and B is the phase constant. The re- 
duction in magnitude of the wave as it travels down the pipeline may be de- 
termined from the attenuation coefficient: 


(5) 


Beta determines the change in phase between the pressure and volume 
current as the pressure wave travels in the pipeline: 


The reader is referred to the original work for the complete derivation of 
the above equations. The primary interest herein is the application of this 
method of analysis to a full-scale field installation. 

Eqs. 1 and 2 will be recognized as being of the same form as the classi- 
cal wave equations of acoustics and electrical engineering, and an electrical- 
hydraulic analogy becomes evident. Therefore: 


Py = & & (7) 


i.e., ohms law for the hydraulic circuit becomes the defining equation for the 
impedance (Z) to the alternating flow. Eqs. 1 through 7 from the basis of the 
analysis of a pumping system. The recognition of the electrical analogy 
permits the formulation of a method of hydraulic circuit analysis similar in 
many respects to electric circuit analysis. 


The Pump Equation 


Inspection of Eq. 7 reveals that, once any two of the three quantities are 
known at a particular point in the pipeline system, the third may be calcu- 
lated. Then, by application of Eqs. 1 and 2, the fluctuating pressure and al- 
ternating volume current may be calculated at other points in the system. 
Pressure is a measurable quantity. The volume current and impedance are 
not. 

However, the use of a reciprocating pump on a pipeline system gives rise 
to an alternating periodic volume current or volumetric flow rate. Since the 
flow is periodic, it may be analyzed by Fourier Series methods and ex- 
pressed in terms of a trigonometric series involving harmonics of a basic 
frequency. The response of the system to the various harmonics may then 
be determined; and, using superposition, the effects of all the harmonics may 
be found if desired. 

When the connecting rod length, ¢, is not relatively long compared to the 
crank arm, r, the piston motion is not simple harmonic, and the piston 
velocity for constant w becomes: 
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= = rw [sine + Sin | (8) 
; The variable flow rate is then: e 
2 
: where A¢ is the area of the cylinder. For a double-acting cylinder A. = 7 - 
n(D? - d2) 
4 for 0<wt < 180° and A, = for 180°<w t < 360° where “D” is the 
-. liner diameter and “d” is the diameter of the piston rod. Figure No. 1 isa ig 
4 plot of the discharge from one double-acting cylinder. ‘% 


T 


4 Fig. 1. Discharge from One Double-Acting Cylinder. * 
a The Fourier Series solution for this case, referred to head end dead x 
a center, is: 
a + Sin 2wt - Cos ewe] 
4 Cos mWt 
m- - 4 
(10) 
bs The effect of referring Eq. 10 to crank end dead center is to reverse the i 
‘ algebraic sign of the odd harmonics. i 
ps For two double-acting cylinders, i.e., a double-acting duplex pump such as a 


| 

; 0 = 

4 
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the Gaso pumps at the Test Station with one cylinder ninety degrees ahead of 
the other, as shown in Figure No. 2, and with the first cylinder referred to 


Fig. 2. Discharge from Double-Acting Duplex Pump. 


head end dead center, the Fourier Series analysis gives: 


Q = | Sin(wt + 135°) 


°) rae 4 
wt 45 s(n 
+ in( + + Sat  gocy | Co (m wt + 45 ) 


mo - 


Cos(mwt + 135°) - 
J 


9, 13 


(11) 


Similar equations may be developed for any multicylinder double -acting 
pump referenced to head end or crank end dead center by superimposing the 
effects of the required number of double-acting single cylinders (Eq. 10) 
properly phased. 

It should be noted that the first term in Eqs. 10 and 11 is the mean flow (q) 
component, while all the other terms are the variable flow components. The 
magnitude of all odd harmonics above the third and all even harmonics above 
the eighth is normally small in comparison to the predominant fourth har- 
monic. Neglecting these terms, Eq. 11 may be reduced to: 
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q Tew ype - 24° 2 r V2 q, 
sin(wt + 4 °)| 4 ve 
5°) + Cos(3wt + 45>) 
n( + 
4 15 3 
: The following quantities are applicable to the Gaso pumps at the Test 
4 Station: 
— D = 5-1/2 inches = 0.458 ft. 
4 d = 2 inches = 0.167 ft. 
r = 6 inches = 0.5 ft. 
4 W = (rpm)&2 = 0.105 (rpm) 
{ = 33-1/2 inches = 2.742 ft. 
ag 
a Eq. 1la then reduces to: 
4 Q@ = 0.0103 (rpm) 0.000816 (rpm) Sin (wt - 128.9°) 
+ 0.000221 (rpm) Sin (3wt + 135°) 
4 0.00137 (rpm) Sin (4wt - 90°) 
4 0.000325 (rpm) Sin (8wt - 90°) 
J Eq. 11b is used throughout the analyses in the following chapters as the 
q value of the variable volume current, Qg, generated by the pump. 
E For a pump with no return rod in the cylinder (i.e., d = 0), Eq. 10 for one 
a double -acting cylinder becomes: 
q- mw |v, 
4 
: E and Eq. 11, for two double-acting cylinders 90° apart, reduces to: 


Equations may be developed for any multicylinder double-acting pump by 
superimposing the effects of the required number of double-acting cylinders 
(Eq. 12) properly phased. 
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Field Tests 


The Test Station was equipped with four double-acting duplex, Gaso pumps 
with V-belt drive. Crude oil is pumped from tankage through approximately 
forty-five miles of eight-inch line into tankage. There are several diameter 
changes and short stub lines involved. 


Instrumentation 


Fluctuating pressures were measured and recorded with the aid of Con- 
solidated Engineering Corporation and Baldwin pressure cells. The output of 
the cells was fed into a Midwestern recording oscillograph. The measured 
values of pressure recorded in Tables II through VII were obtained from the 
analysis of the oscillograph traces. Brown recording pressure gages were 
used to measure the mean discharge pressures at the sending and receiving 
ends. 


Crude Oil Data 


Samples of crude oil were taken before, during, and after the test runs. 
Analysis of the samples gave the following data: 


Density 40.7° API at 54° F. or 1.59 slugs per cubic foot. 
Viscosity 51.5 SUS or 13.45 x 10° Ib. sec. per square foot. 


Velocity of Propagation 


The velocity of propagation of the pressure waves in the system was de- 
termined by measuring the time required for a down surge to travel from the 
sending end pumps to the receiving end. An average time of travel of 59.3 
seconds was used, giving a velocity of propagation of 4029 feet per second. 


Determination of the Attenuation Coefficient 


The attenuation coefficient may be readily determined from Eq. 5. Per- 
tinent data are tabulated in Table I. The mean pressure drop was plotted 
against the mean velocity and a value of 1.75 determined for “n.” The re- 
sulting values of the attenuation coefficient were small enough to be ignored 
in the analysis of the Test Station piping. 


Analysis 


With the distance, x, measured from the receiving end, the pressure and 
volume current for any end termination, Z,;, are given by Eqs. 1 and 2. 
With ax small (i.e., Coshax = 1 and sinhax = 0), Eqs. 1 and 2 reduce to: 


Py, = PypCos@x + (14) 
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& = ZyCosA@x + Sin@x (15) 
Phe 
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At the receiving end: 
(16) 
TABLE I 
DATA FOR DETERMINATION OF n AND 
ax 1076 
Numeric Np 
per ft. x 10° 
111.0 .800 
138.0 129.7 289. 945 
67.0 59.5 
60.0 52.7 
90.0 82.5 


meen pressure as shown by recorder at the sending end. 
mean pressure as shown by recorder at the receiving end. 


Pp PS less PR plus effects of elevation change between stations. 


Letting x = ¢ , the sending end impedance becomes: 


Zo + j2pTanZd Qs 


With the distance, x, measured from the sending end, the pressure and 
volume current for any end termination, Zs, are given as: 


Px = P,CoshYx -. Q,Z,Sinh¥ x (18) 


& = QCosh¥x - Sinn ¥ x (19) 


Again, when ax is small (i.e., Coshax = 1 and Sinhax = 0), Eqs. 18 and 19 re- 
duce to: 


Py P,Cos@x - jQ,Z,Sin@x (20) 


P 
- iz SinB x (21) 
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Letting x = ( and 
(22) 


the receiving end impedance, Z,, becomes: 


Z, - (23) 


- 


The field test pipeline system was basically composed of pumps, pipeline, 
and receiving end tank, as shown in Figure No. 3a. 


Sending End Receiving End 
Pumps Tanks 
Pipeline 


Fig. 3. Analysis Sketches. 


Electrically, this system may be diagramed as shown in Figure No. 3b. 
Analysis of the system based on this arrangement gave pressure magnitudes 
in good agreement with measured values. However, the resulting phase 
angles were generally not in good agreement with measured values. 

The concept of treating the pump as a constant current generator (electri- 
cally speaking) was then introduced. That is, the pump was thought of as 
having an internal impedance, Zj, in parallel and generating a volume current, 
Q,, which divides—part, Q;, circulating through the internal impedance and 
the remainder, Qg, flowing out into the pipeline system proper. See Figure 
No. 3c. 

Then: 


P, = Z,Q,=2,Q, - 2,9, (24) 
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Q, = & + Q, (25) 


From Eqs. 24 and 25: " 


Q Z; Q; 9 

Zp 1i ii (26) 

Qe Qi + + Zs, 

or 

2.2. 

(27) 

* Zp 


The sending end impedance, Z,, may be readily calculated from Eq. 17. 
The measurement of Ps permits calculation of Q, from Eq. 24. Knowing Q, 
and Q, (from Pump Equation), Qj may be easily determined from Eq. 25. 
With the values of Qj and Q, known, the values of Zj and Zp) may be estab- 
lished from Eq. 24 or Eqs. 26 and 27. 

If P, and Z, are known, the value of Q, may be calculated. It is not neces- 
sary to determine the value of Zj in order to proceed with the analysis. 

However, if P, is now known, it becomes necessary to evaluate Z; by some 
means. This, in turn, permits evaluation of Z, (impedance at the pump) by 
use of Eq. 26. Then, the sending end pressure value, Ps, may be calculated 
from Eq. 24. 

In the case of the field tests, Ps was a known (measured) quantity. It was 
then possible to calculate Q, and, hence, values of pressure and volume cur- 
rent at other points in the system. 


One-Pump Tests with Pump at the End of the Line 


The physical arrangement of the equipment was as shown in Figure No. 4. 
Points P2, P3, etc., denote locations of pressure pickups. 


Pump No. 4 


Pe P5 


To Tanks (45 Miles) ——» ) 


Figure No. 4 - Pump at End of Line 


Three tests (Nos. 1, 2, and 3) at 66.7, 77.4, and 82.2 r.p.m. were run using 
this setup. The circuit may be visualized electrically as shown in Figure 
No. 3. 

Applying the equations and techniques outlined in the previous section, the 
pressures at the various pickup points were calculated. Measured and calcu- 
lated values are compared in Tables II, I, and IV. 
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TABLE II: Test No. 1” 


COMPARISON OF MEASURED AND CALCULATED PRESSURE VALUES 


(at Pump) 
Measured 


Harmonic psi Q deg. 


(at Pump) 
Calculated 


c.f.s. Q deg. 
- 90. 
+145. 


- 90. 


Pe (x = 111 feet) 


Calculated Measured 
Harmonic psi Q deg. psi Q deg. 


12.76 -105. 


P3 (x = 32 feet) 
Calculated Measured 


psi Q deg. psi Q deg. 


12.77 - 96.% 2.50 -100.5 


6.2 +127. +122.2 
-118.6 


-150.4 


P5 (x = 554 feet) 


Calculated Measured 


psi Q deg. psi Q deg. 


-159.7 -161.1 


-115.6 


a 
x 
oat 
3 Oe 46 
8 5.16 -104.1 0.0304 -130.5 4.26 5.17 
af. 
Prms. 15.79 15.80 14.90 
3 6.23 + 99.8 4.97 + 90.0 6.23 4.95 - 92.6 
* 
4 16.25 -148.2 15.61 -156.8 16-57 = 12-27. = BIG 
i 8 4.40 4.80 +143.8 4.25 + 2.6 4.12 + 33.7 
Prms. 15.58 14.54 15.13 12.04 
i Pump No. 4 @ 77.4 r.p.m. ae . 
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* 
a TABLE III: Test No. 2 4 
COMPARISON OF MEASURED AND CALCULATED PRESSURE VALUES 
j ; P4 (at Pump) Qs (at Pump) P3 (x = 32 feet 4 
a Measured Calculated Calculated Measured 
= Harmonic psi g deg. c.f.s. Q deg. psi Q deg. psi 9 deg. 
1 12.73 - 94.7 0.0988 - 99.1 12.67 -101.4 12.39 - 98.8 aig 
3 7-17 +140.4 0.0529 +127.5 6.70 +125.2 6.55 +117.2 
15.86 - 64.3 0.1129 -101.0 14.82 -110.2 15.17 -118.6 
8 5.15 -103.4 0.0293 -130.5 4.15 --153-4 5-72 -155-7 
Prms. 15.68 14.87 15.15 
4 Pe (x = 111 feet) P5 (x = 554 feet) hs 
2 Calculated Measured Calculated Measured Bic) 
Harmonic psi deg. psi deg. psi deg. psi deg. an 
4 
1 12.68 12.66 -108.7 12.67 -169.2 10.97 -169.9 
3 6.70 + 96.2 6.82 + 91.3 6.70 78.2 5.94  -102.8 
ve 4 14.82 -148.9 13.12 -151.6 14.67 - 21-4 14.46 - 42.9 
3 3 4.15 4129.2 5.09 +128.6 3.71 - Tob 
Prms. 14.87 14.23 14.79 13.76 
F * Pump No. 4 @ 82.2 r.p.m. 
-16- 
ad 


4 
| 


Pk (at Pump) Q, (at Pump) 
Measured Calculated 
psi deg. c.f.s. deg. 
12.22 -115.9 0.0950 -119.8 
5.95 +153.5 0.0448 +142.8 
13.13 - 83.4 0.0963 - 97.4 
4.95 - 79-1 0.0310 -103.2 


13.82 


Pe (x = 1ll feet) 


Calculated Measured 
12.17 -129.2 13.76 -124.2 
5.80 +113.7 5.84 4111.5 
12.53 -136.1 15.32 -136.7 
4.23  +178.6 4.39 +169.7 
13.35 15.45 


-17- 
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TABLE IV: Test No. 3° 


COMPARISON OF MEASURED AND CALCULATED PRESSURE VALUES 


P3 (x = 32 feet) 


Calculated Measured 
psi Bdeg. fei deg. 
12.17 -121.4 13.05 -118.8 
5.80 +137.2 5.50 +143.6 
12.52 -104.7 13.55 -106.2 
4.23 -118.6 4.50 -120.4 


13.35 14.22 


P5 (x = 554 feet) 


Calculated Measured 


psi Q deg. psi deg. 


12.17 -176.4 12.17 -151.2 
5.80 -27.8 5.71 - 50.5 
12.52 + 35.3 12.40 + &.8 
4.23 +161.4 4.00 +133.5 
13.35 13.24 
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One-Pump Tests with Pump Not at the End of the Line 


The physical arrangement of the equipment was as shown in Figure No. 5. 


ump N 


2m 


Figure No. 5 - Pump Not at End of Line 


Three tests (Nos. 4, 5, and 6) at 59.4, 76.4, and 84.5 r.p.m. were run using 
this setup. The circuit may be visualized electrically as shown in Figure 
No. 3. However, Zg in this case is the parallel combination of the long line 
and the short stub in series with the short six-inch line. 

Again, applying the equations and techniques previously outlined, the pres- 
sures at the various pickup points were calculated. Measured and calculated 
pressure values are compared in Tables V, VI, and VIL. 


Analysis of Test Results 


There were six one-pump tests. The oscillograph traces were analyzed 
for the first, third, fourth, and eighth harmonics at three different points in 
the pipeline. 

In order for the theory to predict accurately the experimental pressure, 
the difference “Y” between the observed pressure and theoretical pressure, 
when analyzed statistically, must have reasonable values when compared to 
the precision of measurement. The differences Y = Pth - Pexp are recorded 
in Table VIL 

When these 72 values are analyzed, it is found that: 


1. “Y” is normally distributed with mean, Y = -0.15 psi = 0. 


=\271 
= 
2. The standard deviation = [a | 2 =1.33 psi. 
3. Therefore, with 95 per cent confidence (i.e., 19-1 odds), the experiment- 
al pressure amplitudes are predicted by the theory within 2.66 psi (two 
standard deviations) of the theoretical pressure, without regard to the 
precision of measurement. It is noted that only one measurement dif- 


fered by more than +4 psi and only eight by more than ¢2 psi. (See 
Table VII.) 
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TABLE V: Test No. 4” 


COMPARISON OF MEASURED AND CALCULATED PRESSURE VALUES 


Pl (at Pump) Qs (at Pump) Pe (x = 32 feet) 
Measured Calculated Calculated Measured 


deg. c.f.8. Q deg. psi Q deg. psi Q deg. 


0.0687 8.72 - 94.9 8.61 
0.0378 126. , +123.4 

0.0920 -129.8 10.23 
0.0312 8. -160.9 1.79 


9-99 


P3 (x = 111 feet) P5 (x = 507 feet) 


Calculated Measured Calculated Measured 
Harmonic psi deg. psi deg. psi Bdeg- psi deg. 
8.80 8.94 - 95.2 8.68 -136.4 7-89 -140.4 
5-70 23 4.62 +101.8 5-30 - 2.6 3.19 - 49.4 
11.41 2 12.17 -128.7 9.94 ; 8.86 + 54.5 
3.90 -159.9 2.07 1.86 -162.1 


11.30 10.16 8.78 


* Pump No. 1 @ 59.4 r.p.m. 


if 
hy 
af 
é 
3 4.18 +117.5 + 99.6 
4 4 8.54  -106.5 -130.6 
8 1.40 - 56.3 -166.9 
ad 
-19- 
as 
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Harmonic psi 
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Prms. 


Harmonic psi 
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Prms. 
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TABLE VI: Test No. 5” 
COMPARISON OF MEASURED AND CALCULATED PRESSURE VALUES 
Pl (at Pump) Q, (at Pump) P2 (x = 32 feet 
Measured Calculated Calculated Measured 

Q deg. c.f.s. deg. psi deg psi deg. 
11.92 -131.2 0.0962 -126.2 12.09 -137.6 12.48 -141.0 
3-77 +136.2 0.0395 +146.6 4.320 4124.1 4.39 74.8 
7-21 -108.4 0.0969 -100.2 9.27 -143.7 12.83 -161.3 
4.10 -21.4 0.0472 -106.0 1.59 +176.3 1.14 +146.0 

10.60 11.26 13.05 

P3 (x = lil feet) P x = 507 feet 

Calculated Measured Calculated Measured 

Q deg. psi deg. psi deg. psi deg. 
12.25 +-137-6 13.72 -133.5 12.09 -191.6 12.7%  -174.2 
-+114.1 5.28 + 87.2 4.32 47.9 3-77 -101.6 
11.73 -143.7 16.67 -148 .6 S27 9.77 <- 
6.48 +176.3 6.28 +163.3 1.59 +104.1 5.09 + 85.5 

13.30 16.33 11.25 11.97 


* Pump No. 1 @ 76.4 r.p.m. 
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TABLE VII: Test No. 6” 


COMPARISON OF MEASURED AND CALCULATED PRESSURE VALUES 


Pl (at Pump) 
Measured 


psi deg. 


Q, (at Pump) 
Calculated 


c.f.s. Q deg. 


Pe (x = 32 feet) 


Calculated Measured 


psi Q deg. psi Q deg. 


Harmonic 


-114.9 0.0639 -109.4 8.01 -122.0 8.21 -125.9 


4.62 + 97-2 


-142.7 


+115.7 0.0437 +126.3 + 90.2 5.90 


-104.3 0.1159 - 99.8 10.59 -146.4 12.14 


- 42.6 0.0629 -131.6 0.83 +146.4 0.23  -135.0 


9.96 11.20 


P x = 1ll feet P x = 507 feet 
Calculated Measured Calculated Measured 


psi Q deg. psi Q deg. psi y deg. psi . Q deg. 


Harmonic 


8.15 
5.40 
14.18 
7.38 
13.25 


* 
Pump No. 1 @ 84.5 r.p.m. 


-122.0 
+ 90.2 
-146.4 
+146.4 


9.61 
6.78 
15.39 
7-94 
14.80 


-131.7 
+ 94.8 
-146.2 


+152.9 


8.01 
4.61 
10.59 
0.83 
9.96 


-181.2 
89.5 
25.4 
28.4 


10.40 
4.69 
9.94 
0.51 


10.71 


-178.1 
- 99-9 
- 20.4 
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TABLE VIII 


DIFFERENCE BETWEEN MEASURED AND CALCULATED 


PRESSURE VALUES FOR EACH HARMONIC 


Y 
Difference 


Pickup Harmonic psi Test Pickup 


+ 


OFWHKOFWH 
+ 


' 

' 


ae) 
\O 


Ww 
Wr OFW 


or 
Ww oO 
Www 


Plus indicates calculated value greatest. 


-23- 


Harmonic 


OFW 


Y 
Difference 


psi 
DSi 


‘ 
q 
| March, 1957 
l +0.29 
+0.15 
-0.35 
-1.57 
2 +0.02 a 
-0.12 
+1.70 
-0.94 
5 +1.70 
+0.76 
+0.21 
+0.44 
3 3 2 +0.11 
| +1.08 
-0.29 af 
+0.24 
| 3 -0.14 
| +1 -08 
-0.76 
+0.13 
+2.11 
+1.00 
+0.21 
+0.20 
+0.60 
-1.46 
-0.56 
-2.39 
-0.08 
+0.65 
+0.32 
i 
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TABLE IX 


MPARISON OF MEASURED AND CALCULATED 
ROOT MEAN SQUARE PRESSURE VALUES 


Calculated Measured 
Pickup Prms Prms Difference 


No. psi psi* 


14.90 +0.90 
14.54 +1.04 
12.04 +3.09 


Ww Ww 


15.15 -0.28 
14.23 +0.64 
13.76 +1.03 


fo 


14. -0.87 
15. -2. 
+0. 


9. +0. 
11.4 -O. 
8. +1. 


3 
3 
3 


ww fo 


13. -l. 
16. -3. 
11. -0. 


Vw 


11.2 
14.8 
10. -0. 


ON ON ODN 


* Plus indicates calculated value greater. 
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When considering that the factors that comprise the error of measurement 
alone might be as much as +4.32 psi, one can only conclude the fundamental 
principles and theory satisfactorily predict the measured pressures well 
within the limits of the instrumentation. 

In comparing the root mean square (square root of one-half the sum of the 
Squares) pressure values, it was found that none of the eighteen values ex- 
ceeded the probable error of +4.32 psi. (See Table IX.) 


CONCLUSIONS AND RECOMMENDATIONS 


The results of the work related herein adequately account for the pressure 
fluctuations in a reciprocating pump system. Specifically: 


1. The validity of the basic underlying derivations has been established. 
a) The fundamental principles and equations formulated in the basic 
theoretical analysis have been verified. 
b) The principle of superposition is applicable to this situation. 
c) Simplification of the general mathematical expressions is valid for 
this case. 


The applicability of the fundamental principles to a field-scale system 

has been demonstrated. 

a) The techniques of analysis were successfully applied to an actual 
field installation. 

b) It appears that the use of these results will be limited only by the 
ingenuity of the persons performing the analysis. 

c) The results of this work are sufficiently general enough in nature to 

permit the designer to apply them to a variety of piping configura- 

tions with a high degree of confidence. 


The parameters based on the concepts of the hydraulic-electric analogy 
have been reasonably well established. Although work should be continued on 
the study of these basic parameters, the main effort should be directed to- 
ward further development of techniques of analysis for more complex pipeline 
systems. 

The basic theme of this entire research activity is the hydraulic-electric 
analogy. This being the case, electronic computer techniques could be ap- 
plied and would greatly facilitate the preparation of much of the material in 
chart and graph form. Such procedures could lessen considerably the com- 
putational effort required in the analysis of many pressure surge problems. 
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APPENDIX - NOTATION 


The symbols listed below are briefly defined, and the descriptive units 
refer to the gravitational (technical) system. “F” is a force usually in 
pounds, “L” is a length usually feet, and “T” is time usually seconds. The 
mass unit is the slug = 32.2 lb/g. 


- The velocity of propagation of the pressure wave. LT71, or ft. 
per sec. 

- Area. L2, square ft. 

As a subscript refers to a particular point in pipeline system; 

see appropriate sketches. 

- Asa subscript refers to the stub or branch in the pipeline 
system. 

- The diameter of the piston rod. L, ft. 

- The diameter of the pipe for the pipeline analysis. L, ft. 

The diameter of the liner for the pump analysis. L, ft. 

- Asa subscript refers to equivalent parallel combination of 
branch and main line impedances. 

- The frequency. T-1, cycles per second. 

- The gravitation constant. LT-2, 32.2 ft. per sec2, 

- Asa subscript refers to internal mechanism of pump. 

- 1, i.e., j is complex operator = V-1. 

Length. L, ft. 

- An integer. 

- Asa subscript refers to the main line. 

- The exponent of U for mean flow. 

- The Reynold’s Number. 

- Asa subscript refers to the equivalent impedance at the pump 
(i.e., parallel combination of sending end and internal pump im- 
pedances). 

p - The — yn | causes a mean flow in a finite length of 
line. lb. per ft.2, 

Psub - The anuatin representation of the standing wave pressure am- 
a. The subscript denotes the point in question. FL~2, lb. 
per ft.¢. 

q - The mean rate of flow. L°T- 1, cubic ft. per sec. 

Qsub - The vectorial representation of the volume current or harmonic 
flow rate. L3T-1, cubic ft. per sec. 

Qg - The volume current generated by the oscillating piston of the re- 
ciprocating pump. L3T-1, cubic ft. per sec. 

r - Asa subscript refers to the receiving end of the pipeline. 

r - The radius of the crank arm or the eccentric of the crankshaft. 
L, ft. or 1/2 of stroke. 

s - As a subscript refers to the sending end of the pipeline. This is 
usually at the pump end. 

s - The stroke of the piston (2r). L, ft. 
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Fa - Time. T, sec. 
U - The mean velocity of the fluid in the pipeline. LT-1, ft. per sec. 
Yp - The velocity of the piston in a reciprocating pump. LT-|, ft. per 
sec. 
x - The distance from the receiving end or sending end. See appro- 
priate sketches. L, ft. 
x - Used as a second subscript to denote existence of several 
harmonics. 
Ze - The characteristic impedance of the pipeline. FL-5T, lb. sec. a 
per ft.5, 
Zsub - The complex quotient of the hydromotive pressure (P) and hy- om 


draulic volume current (Q) at any point in the system denoted by 
the subscript. FL-5T, lb. sec. per ft.9. 

a - The attenuation constant. L-1, numeric per ft. 

B - The phase constant. L-1, radians per ft. or degrees per ft. 

Y - The propagation constanta + B . L-1, 
p 


- The density of the fluid. FL-1T2 per L3 or FL-4T2, slugs per 
ft. 

- The phase angle of the polar form of the vectors (degrees). 

- The circular speed of the rotating vectors; also the angular 
velocity of the crankshaft. L-1, radians per second. 

This sign as used herein means approximately equals. 


NOTE: All capital P’s, Q’s, and Z’s in equations in this report refer to 
vector quantities. When a single subscript is used, it refers toa 

particular position. When two subscripts are used, the first refers 
to position, the second to the particular harmonic in question. 
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PIPELINE DIVISION 


Proceedings of the American Society of Civil Engineers 


NEWS 
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The PIPELINE DIVISION was inaugurated at the annual convention in 
Pittsburgh last October. At these first meetings each chairman of the ad- 
ministrative and technical committees was assigned the responsibility for 
organizing his committee to work with Local Sections, the Junior Members 
and Student Chapters. Also, to correlate his task groups with work done by 
other divisions of the Society, the API, AGA, and ASME in particular, other 
societies and organizations, to avoid any duplication of effort. At the 
Jackson, Mississippi, convention in February, task groups with schedules for 
this year’s accomplishments were reported by all chairmen. 

Local Sections have been contacted requesting their nominations of one or 
more co-ordinators for our Pipeline Division activities. Also, for a list of 
members who desire affiliation with our division, since this information is 
not available at Headquarters. This will add to our mailing list for our 
JOURNAL and division notices. If you know of anyone who should affiliate 
with us, please urge them to advise Headquarters direct. 

Volunteer workers are still needed for all of the administrative and 
technical committees. If you are interested, please communicate with the 
committee chairmen listed below, or any member of the Executive Committee: 


Publications - H. A. Smallwood, Assoc. Professor of Civil Engineering, 
Madison, Wisconsin. 


Session Programs - R. E. Kling, Chief Engineer, Site Oil Company, 3010 
Locust Street, St. Louis 3, Missouri. 


Membership - L. M. Odom, Pyburn & Odom, Consultants, P. O. Box 267, 
Baton Rouge, Louisiana. 


Public Relations - D. M. Taylor, Gulf Coast Editor, “The Petroleum Engi- 
neer,” Rice Professional Bldg., 2370 Rice Blvd., Houston, 
Texas. 


Cooperation With Local Sections - R. R. Crawford, Senior Engineer, Bechtel 
Corp., 209 H Street, San Rafael, California. 


Fluid Dynamics - W. T. Ivey, Statistician & Planning Engineer, Southern 
Natural Gas Company, P. O. Box 2563, Birmingham 2, Alabama. 


Note: No. 1957-6 is part of the copyrighted Journal of the Pipelines Division of the 
American Society of Civil Engineers, Vol. 83, PL 1, March, 1957. 
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Pipeline Crossings of Railroads & Highways - J. E. Thompson, Supt. of Pipe- 
line Construction, Natural Gas Pipeline Co. of America, 122 
South Michigan Ave., Chicago 3, Dlinois. 


Pipeline Location, Surveying & Mapping - H. M. Hayes, Profect Engineer, 
Ford, Bacon & Davis, Inc., P. O. Box 1762, Monroe, Louisiana. 


Pipeline Design, Specifications & Operating Standards - C. L. Shea, Senior 
Project Engineer, Products Pipelines, Shell Oil Company, 
4459 North Keystone Avenue, Indianapolis 5, Indiana. 


Pumping & Compressor Stations - F. E. Culvern, Research Engineer, 
Panhandle Eastern Pipeline Co., 1221 Baltimore Ave., Kansas 
City 5, Mo. 


Storage of Pipeline Fluids - chairman to be appointed. 
Executive Committee: 


Chairman E. V. Hunt, Chief Engineer, The Alberta Gas Trunkline Co. 
Ltd., 320 - 9th Avenue West, Calgary, Alberta, Canada. 

Vice Chairman A. E. Poole, President, The Hallen Construction Co., 
Inc., 4270 Austin Boulevard, Island Park, New York. 

S. E. Huey, S. E. Huey & Company, Heninger Building, Monroe, 
Louisiana. 

F. C. Culpepper, Jr., Project Manager, Ford, Bacon & Davis Construc- 
tion Corporation, P. O. Box 710, Alliance, Ohio. 

Secretary J. B. Spangler, Transcontinental Gas Pipe Line Corporation, 
P. O. Box 296, Houston 1, Texas. 


Contact Member from Board of Direction, L. A. Elsener, V. P., Chicago 
Bridge & Iron Company, 100 Bush Street, San Francisco 4, 
Calif. 
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PROCEEDINGS PAPERS 


The technical papers published in the past year are identified by number below. Technical - 
division sponsorship is indicated by an abbreviation at the end of each Paper Number, the 
symbols referring to: Air Transport (AT), City Planning (CP), Construction (CO), Engineering 
Mechanics (EM), Highway (HW), Hydraulics (HY), Irrigation and Drainage (IR), Power (PO), 
Sanitary Engineering (SA), Soil Mechanics and Foundations (SM), Structural (ST), Surveying and 
Mapping (SU), and Waterways and Harbors (WW) divisions. Papers sponsored by the Board of 
Direction are identified by the symbols (BD). For titles and order coupons, refer to the appro- 
priate issue of “Civil Engineering.” Beginning with Volume 82 (January 1956) papers were 
published in Journals of the various Technical Divisions. To locate papers in the Journals, the 
symbols after the paper numbers are followed bya numeral designating the issue of a particular 
Journal in which the paper appeared. For example, Paper 1113 is identified as 1113 (HY6) 
which indicates that the paper is contained in the sixth issue of the Journal of the Hydraulics 


Division during 1956. 
VOLUME 82 (1956) 


MARCH: 906(WW1), 907(WW1), 908(WW1), 909(WW1), 910(WW1), 911(WW1), 912(WW1), 913 
(WWI1)©, 914(ST2), 915(ST2), 916(ST2), 917(ST2), 918(ST2), 919(ST2), 920(ST2), 921(SU1), 
922(SU1), 923(SU1), 924(ST2)¢. 


APRIL: 925(WW2), 926(WW2), 927(WW2), 928(SA2), 929(SA2), 930(SA2), 931(SA2), 932(SA2)°, 
933(SM2), 934(SM2), 935(WW2), 936(WW2), 937(WW2), 938(WW2), 939(WW2), 940(SM2), 941 
(SM2), 942(SM2)°, 943(EM2), 944(EM2), 945(EM2), 946(EM2)°, 947(PO2), 948(PO2), 949(PO2), 
950(PO2), 951(PO2), 952(PO2)¢, S53(HY2), 954(HY2), 955(HY2)°, 956(HY2), 957(HY2), 958 
(SA2), 959(PO2), 960(PO2). 


MAY: 961(IR2), 962(IR2), 963(C P2), 964(CP2),965(WW3), 966(WW3), 967(WW3), 968(WW3), 969 
(WW3), 970(ST3), 971(ST3), 972(ST3)°, 973(ST3), 974(ST3), 975(WW3), 976(WW3), 977(IR2), 
978(AT2), 979(AT2), 980(AT2), 981(IR2), 982(IR2)¢ ,983(HW2), 984(HW2), 985(HW2)°, 986(ST3), 
987(AT2), 988(CP2), 989(AT2). 


JUNE: 990(PO3), 991(PO3), 992/PO3), 993(PO3), 994(PO3), 995(PO3), 996(PO3), 997(PO3), 998 
(SA3), 999(SA3), 1000(SA3), 1001(SA3), 1002(SA3), 1003(SA3)©, 1004(HY3), 1005(HY3), 1006 
(HY3), 1007(HY3), 1008 (HY3), 1009 (HY3), 1010 (HY3)°, 1011 (PO3)°, 1012 (SA3), 1013 (SA3), 
1014(SA3), 1015(HY3), 1016(SA3), 1017(PO3), 1018(PO3). 


JULY: 1019(ST4), 1020(ST4), 1021(ST4), 1022(ST4), 1023(ST4), 1024(ST4)©, 1025(SM3), 1026 
(SM3), 1027(SM3), 1028(SM3)°, 1029(EM3), 1030(EM3), 1031(EM3), 1032(EM3), 1033(EM3)°. 


AUGUST: 1034(HY4), 1035(HY4), 1036(HY4), 1037(HY4), 1038(HY4), 1039(@HY4), 1040(HY4), 
1041(HY4)¢, 1042(PO4), 1043(PO4), 1044(PO4), 1045(PO4), 1046(PO4)¢, 1047(SA4), 1048 
(SA4)¢, 1049(SA4), 1050(SA4), 1051(SA4), 1052(HY4), 1053(SA4). 


SEPTEMBER: 1054(ST5), 1055(ST5), 1056(STS), 1057(ST5), 1058(STS), 1059(WW4), 1060(WW4), 
1061(WW4), 1062(WW4), 1063(WW4), 1064(SU2), 1065(SU2), 1066(SU2)°, 1067(sTS)“, 1068 
(wwa4)°, 1069(WWw4). 


OCTOBER: 1070(EM4), 1071(EM4), 1072(EM4), 1073(EM4), 1074(HW3), 1075 (HW3), 1076(HW3), 
1077(HY5), 1078(SA5), 1079(SM4), 1080(SM4), 1081(SM4), 1082(HY5), 1083(SA5), 1084(SA5), 
1085(SA5), 1086(PO5), 1087(SA5), 1088(SA5), 1089(SA5), 1090(HW3), 1091(EM4)°, 1092 
(HYS)©, 1093(HW3)*, 1094(POS5)°, 1095 (sm4)°. 


NOVEMBER: 1096(ST6), 1097(ST6), 1098(ST6), 1099(ST6), LIOL(ST6), 1102(IR3), 
RS), NOS(IR3), 106(ST6), 107(ST6), 1108(ST6), 110(AT3)°, 


DECEMBER: 1113(HY6), 1114(HY6), 1115(SA6), 1116(SA6), 1117(SU3), 1118(SU3), 1119(WW5), 
1120(WW5), 1121(WW5), 1122(WW5), 1123(WW5), 1124(WW5)°, 1125(BD1)“, 1126(SA6), 1127 
(SA6), 1128(WW5), 1129(SA6)©, 1130(PO6)°, 1131(HY6)°, 1132(PO6), 1133(PO6), 1134(PO6), 
1135(BD1). 

VOLUME 83 (1957) 


JANUARY: 1136(CP1), 1137(CP1), 1138(EM1), 1139(EM1), 1140(EM1), 1141(EM1), 1142(SM1), 
1143(SM1), 1144(SM1), 1145(SM1), 1146(ST1), 1147(ST1), 1148(ST1), 1149(ST1), 1150(ST1), 
1151(ST1), 1152(CP1)©, 1153(HW1), 1154(EM1)°, 1155(SM1)°, 1156(ST1)°, 1157(EM1), 1158 
(EM1), 1159(SM1), 1160(SM1), 1161(SM1). 


FEBRUARY: 1162(HY1), 1163(HY1), 1164(HY1), 1165(HY1), 1166(HY1), 1167(HY1), 1168(SA1), 
1169(SA1), 1170(SA1), 1171(SA1), 1172(SA1), 1173(SA1), 1174(SA1), 1175(SA1), 1176(SA1), 
1178(SA1), 1179(SA1), 1180(SA1), 1181(SA1), 1182(PO1), 1183(PO1), 1184(PO1), 


MARCH: 1186(ST), 1187(ST), 1188(ST), 1189(ST), 1190(ST), 1191(ST), 1192(ST)°, 1193(PL), 
1194(PL), 1195(PL), 


c. Discussion of several papers, grouped by Divisions. 
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